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ABSTRACT 

Context. Bar-induced gas inflows towards galaxy centres are recognised as a key agent for the secular evolution of galaxies. One 
immediate consequence of this inflow is the accumulation of gas in the centre of galaxies where it can form stars and alter the 
chemical and physical properties. 

Aims. Our aim is to study whether the properties of the ionised gas in the central parts of barred galaxies are altered by the presence 
of a bar and whether the change in central properties is related to bar and/or parent galaxy properties. 

Methods. We use a sample of nearby face-on disc galaxies with available SDSS spectra, morphological decomposition, and informa¬ 
tion on the stellar population of their bulges, to measure the internal Baimer extinction from the Her to Hj3 line ratio, star formation 
rate, and relevant line ratios to diagnose chemical abundances and gas density. 

Results. The distributions of all the parameters analysed (internal Baimer extinction at Up (c(Hp)), star formation rate per unit area, 
electron density, [N ii]/16583/Hcr emission-line ratio, ionisation parameter, and nitrogen-to-oxygen (N/0) abundance ratio) are differ¬ 
ent for barred and unbarred galaxies, except for the R 23 metallicity tracer and the oxygen abundance obtained from photoionisation 
models. The median values of the distributions of these parameters point towards (marginally) larger dust content, star formation rate 
per unit area, electron density, and ionisation parameter in the centres of barred galaxies than in their unbarred counterparts. The 
most remarkable difference between barred and unbarred galaxies appears in the [N ii]/16583/Hcr line ratio that is, on average, ~25% 
higher in barred galaxies, due to an increased N/0 abundance ratio in the centres of these galaxies compared to the unbarred ones. 
We analyse these differences as a function of galaxy morphological type (as traced by bulge-to-disc light ratios and bulge mass), 
total stellar mass, and bulge Sersic index. We observe an enhancement of the differences between central gas properties in barred and 
unbarred galaxies in later-type galaxies or galaxies with less massive bulges. However, the bar seems to have a lower impact on the 
central gas properties for galaxies with bulges above ~ 10 ^^ M© or total mass > 10 ^^-^ M©. 

Conclusions. We And observational evidence that the presence of a galactic bar affects the properties of the ionised gas in the central 
parts of disc galaxies (radii < 0.6-2.1 kpc). The most striking effect is an enhancement in the N/0 abundance ratio. This can be 
interpreted qualitatively in terms of our current knowledge of bar formation and evolution, and of chemical evolution models, as 
being the result of a different star formation history in the centres of barred galaxies caused by the gas inflow induced by the bar. Our 
results lend support to the scenario in which less massive and more massive bulges have different origins or evolutionary processes, 
with the gaseous phase of the former currently having a closer relation to the bars. 

Key words. Galaxies: evolution - (Galaxies:) bulges - Galaxies: spiral - Galaxies: abundances - Galaxies: ISM - ISM: general 


1. Introduction 

The distribution of gas and stars in disc galaxies has a basic ax- 
isymmetric structure, which is driven by rotation. However, su¬ 
perimposed non-axisymmetric morphological features, such as 
spiral arms or bars, may have important consequences on their 
evolution (e.g. jKormendy & Kennicutt|2004| ). This secular, inter¬ 
nally driven, evolution is thought to be dominant in recent times, 
as opposed to the evolution driven by galactic hierarchical merg¬ 
ing and external gas accretion (e.g. [Barnes & Hernqui^|1992| ). 
Stellar bars have received a lot of attention and are seen as im¬ 
portant agents for the secular evolution of galaxies, because their 
non-axisymmetric gravitational potential redistributes the angu- 
lar momentum of the gas and stars in galactic discs (e.g. jSell-j 
wood & Wilkinson|1993[|Athanassoula|2003| ). Gas is dissipative 
and, at radii well inside corotation, loses angular momentum as 


a result of the gravitational torque exerted by the bar. As a result, 
gas is driven towards the galaxy centre (e.g . j Athanassoula| 1 992] 


Friedli & Benz 


|1993t [Friedli et~ar |1994| ). Analytical and nu 


merical simulations also show that the effectiveness with which 
a bar influe nces galaxy dynamic s (bar strength) changes as the 
bar evolves ( [Athanassoula|2003'] ) and, as a consequence, the bar- 
induced gas inflow rate changes accordingly ([Regan & Teuben 
|2004l ). 

According to simulations, bar-induced gas inflow can affect 
the central properties of galaxies in a number of ways: it can 
increase the central gas and dust concentration and trigger star 
formation (SF) ( [Friedli & Benz|199'3l [Martinet & Friedli|1997[), 
produce the formation and growth of a disky pseudobulge ( [Ko- 
[rmendy & Kennicutt[[2004[ [Cheung et al. [2013[ ), and alter the 
central abundance and dilute the initial disc metallicity gradi- 
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ent ( [Friedli & Benz|p^93[ [Friedli et~aL]|1994| ). If the gas were 
to reach the central parsecs it could potentially act as fuel for a 
central massive black hole ( [Shlosman et al.|1989|) and thus form 


an active galact ic nucleus or AGN (e.g. |Coelho & Gadotti|2011 
Oh et al.|2012| ). Eventually, bars can produce central mass con¬ 


centrations and, once they have accumulated sufficient mass in 
the centre, they weaken. For a sufficiently high mass, the bar 
may even d issolve (e.g. [Pfenniger & Norman 199Qt |Friedli &] 
Benz||1993| [Bournaud & Combes ||2002| ), although, in more re¬ 


cent simulations, bars appear more robust than initially though t 
(e.g. |Shen & Sellwood| 2004 [Athanassoula et al.]|2QQ5| |2Q13| ). 


For a detailed des cription of the role of bars in se cular evolution, 
see the review by |Kormendy & Kennicutt| ( [2004| ) and references 
therein. 

Observational evidence exists for bar-induced gas inflows 
(e.g. [Regan et al.|1997t jZurita et ar]|2004| ). The inf low rate has 
been measured to within a range of ~0.1-1 M© yr“^ |Regan et"^ 


1997[ [Sakamoto et al.||1999|), with a higher r ate being found in 


earlier-type galaxies ( [Regan & Teuben|2004t [Sheth et al.|2005[ ) . 

However, the observational search for the theoretically predicted 
effects of bar inflows on the properties of the central gas-phase 
component of galaxies has not been successful in many cases, 
and conclusions are frequently contradictory. 

Compared to unbarred galaxies, barred galaxies have higher 
central molecular mass concentrations ( [Sakamoto et al.[p~999[ 
Sheth et aL][2005| ), which produce a higher star formation rate 


(SFR) in the centres of barred galaxies, as observed by several 
authors (e.g. [Hummel et al. |1990t|Ho et al.[[1997bt[Wang et al.| 


[MT 2 I [Zhou et al.[[2Q15[ [Ellison et al.[[2Ql 1[ [Oh et al.[[2Q12[ ). 
However, there is no agreement on the dependence of this en- 
hancem ent in SFR with bar p roperties. For example, jWang et aT 


( [2Q12[ ) and Zhou et ah] ( [201 5 1 ) find a positive trend between SFR 


[2001 [[Ellison et al.|2011 ). Meanwhile, [Ellison et al.[ ( [2011[ ) find 

an enhancement in central SFR only for barred galaxies with a 
mass greater than 10^^ M©, and the enhancements observed by 
[Oh et al.j p012[ ) are restricted to the reddest galaxies in their 
sample. Moreover, [Cacho et~ar] ( |2014[ ) have recently found no 
enhancement in SFR in barred galaxies compared with unbarred 


ones. 


Studies of the cen tral gas-phase metal licity yielded even 
more complex results. [Ellison et al.[ ( [201 1[ ) find that there are 
central oxygen abundances ~0.06 dex larger in barred than in 
unbarred ga laxies at a given galactic stellar mass, but [Cacho[ 
et al.[ ( [MT^ find no diffe rence in a similar study to the one per¬ 


formed by Ellison et~ar] pOl l) . Other authors ( [Considere et al. 
[2000[ [Dutil & Roy|1999| ) used a smaller galaxy sample and ob 


The advent of the SDSS spectroscopic data has made it pos¬ 
sible to study the effects of bars in central gas properties using 
large galaxy samples and covering wide ranges in host galaxy 
properties. Although studies in this field have contributed to in¬ 
teresting results, all publications, so far, have focused on central 
star formation. Only in a few cases has oxygen abundance also 
been studied ( [Ellison et al.[[2011| [Cacho et al.|2014[ ). However, 
bar-induced gas inflows might also alter other physical proper¬ 
ties of the central gas. For example, if bars are able to produce 
higher centr al mass concentrations ( [Sakamoto et al.|1999|[Sheth[ 
[et al.pOOS] ), higher concentrations of dust should be expected 
in the centres of barred galaxies, which would produce greater 
extinction. This greater concentration of gas and dust in the cen¬ 
tres of barred galaxies, combined with a larger bar-induced SFR, 
could also lead to higher electr on densities in the centres of 
barred galaxies ( [Ho et al.|1997b[ ). The history of star-formation 
in galaxy centres is likely to be different in barred and unbarred 
galaxies due to the gas supply induced by the bars, which might 


and ellipticity or ellipticity-based parameters for strong bars, and 
[Oh et al.[ ( [2012[ ) find it with bar length, but no dependency is 
found by other authors ( Pompea & Rieke[[199Q[ [Roussel et ah] 


tained lower central oxygen abundances in barred starbursts than 
in normal, unbarred, galaxies. So far no relation has been found 
between bar strength and central oxygen abundance. 

Thus, while simulations predict strong effects in the central 
gas properties of barred galaxies, no convincing observational 
confirmation has yet been found. The same bars that produce 
an enhancement in SFR do not seem to produce the metallic- 
ity alteration that simulations predict ( [Friedli et al.|1994[ [Martel[ 
[et al.[[2013[ ). Bar evolution, with its accompanying bar-strength 
change, linked to the fact that bars can even be destroyed (mean¬ 
ing that currently unbarred galaxies have maybe had a bar) might 
be part of the reason for the current discrepancies in findings. 
Also, central properties might also depend on the availability of 
gas and on host galaxy properties in general (namely total mass 
or Hubble-type), which makes sample selection critical. 
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leave footprints in the N/0 abundance ratio (Molla & Gavilan 
[201Q[ [Mallery et al.[[2007l [Edmunds & Pagel[[1978[ ). Conse- 
quently, it is worth exploring other properties of the ionised gas 
that might be modified by the effect of a bar. This is the motiva¬ 
tion for the work presented below. 

In this paper we study ionised gas properties in the centres 
of barred and unbarred galaxies. In addition to the SFR and oxy¬ 
gen abundance, we also measure the Baimer extinction, electron 
density, and the N/0 abundance ratio. We use a sample of nearby 
galaxies with available spectra from which we remove AGN. 
Our galaxy sample has important advantages compared to pre¬ 
vious samples: galaxies are face-on (i < 26°), the stellar com¬ 
ponent has already been studied by [Coelho & Gadotti[ ( [201 1| ), 
and morphological 2D decomposition (and therefore structural 
information on bars, bulges and discs) is also available ( [Gadotti 
[2009[ ). Our aim is to look for observational evidence that bar- 
induced inflows have important consequences for the central 
properties of the gas and, therefore, on secular evolution. 

This paper is organised as follows. In Sect.we describe our 
galaxy sample. In Sect. we describe the procedure to measure 
the emission-line fluxes on which our study is based. Afterwards 
we remove AGN from the sample (Sect.|^. In Sects. Silzl and 
the internal extinction, central oxygen abundance, N/0 ratio, 
SFR, and electron density are calculated for all sample galaxies 
and the distribution of these parameters is analysed and com¬ 
pared separately for both barred and unbarred galaxies. In the 
following sections Em and m we explore the dependence 
of the found differences on Hubble-type, total galaxy, and bulge 
mass, and possible dependencies on bar properties. In Sec t, 
we discuss our results, which are then summarised in Sect.[13[ 


2. Galaxy sample 


The galaxy sample studied here is the one used by [Coelho & 
[Gadotti[ ( [201 1[ ) to study the stellar populations of the bulges of 
barred and unbarred galaxies. This sample is based on the one 
studied morphologically by [GadottT| ( [2009| ). 

The sample contains all spiral face-on galaxies (axial ra¬ 
tio b/a > 0.9) in the Sloan Digital Sky Survey (SDSS) Data 
Release 2, with stellar masses larger than 10^^ M©, redshift 
0.02 < z < 0.07, bulge-to-total luminosity ratio below 0.043 
(i.e. earlier than ~Sd), and with signal-to-noise ratio in their cor¬ 
responding SDSS spectra greater than or equal to 10, measured 
over the spectral range corresponding to the SDSS g-band. This 
sample is made up of 251 barred and 324 unbarred galaxies. The 
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two sub-samples of barred an d unbarred galaxies have similar 
total stellar mass distributions ( [Coelho & Gadotti|2011| ). 


3. Analysis of the spectra 

3.1. Spectra pre-processing 


The spectra analysed in this paper come from the SDSS Data 
Release 7. The spectra pre-processing comprises shifting of the 
original spectra to the rest frame using the redshifts from the 


SDSS database, re-sampling in steps of lA, and Galactic e xtinc 
tion correction as described by |Cid Fernandes et al.| ( |2QQ5] ). 

Afterwards the emission of the stellar compone nt was mod 
elled with the sp ectral synthesis code STARLIGHT ( |Cid Fern^ 


des et al.|2005| ). This is done by comparing the observed SDSS 


spectra, on a pixel-by-pixel basis, to stellar population models 
( [Vazdekis et al. 12010 ). 

We refer the reader to [Coelho & Gadottl] ( |2011| ) and jGadotS] 
( 2009| ) and references therein for full details on the sample selec¬ 
tion and biases, and on the spectra pre-processing and modelling 
of the stellar component. 


3.2. Emission-line measurements 

The emission-line fluxes were automatically measured over the 
stellar-emission-subtracted SDSS spectra of the galaxy sample. 
The latter had already been processed as described in the para¬ 
graph above. 

We developed a code for measuring the fluxes of the 
brightest emission lines of the optical spectra. These 
include [Oii]TT3727,3729, H^, Hy, H/5, [Oiii]T4959, 

[Oiii]T5007, [Oi]T6300, [Nii]T6548, Ha, [Nii]T6584, and 
[S ii]4/16717,6731. All line fluxes were measured independently 
of each other, with no assumptions on the relative strength of 
any pair of lines. 

Initially, the code searches for signiflcant emission (over the 
continuum level) at the rest-frame wavelength of the brightest 
emission lines. At this point, the continuum level is determined 
using a linear fit of the median continuum level in two spectral 
regions on both sides of each emission line. If the code detects 
an emission peak, with a signal of at least twice the standard de¬ 
viation of the continuum level, then the emission line profile is 
fitted with a Gaussian function, and the central wavelength and 
full-width at half maximum (FWHM) of the line are determined. 
Only fitted emission features yielding reasonable value^for the 
central wavelength and FWHM are considered as true detections. 
Afterwards, the code measures the emission-line fluxes. To do 
this, the code first defines new regions dominated by continuum 
emission on both sides of each emission line. The size and dis¬ 
tance of these continuum regions from the emission-line central 
wavelength depend on both the previously determined FWHM 
and on the location of neighbouring emission lines. The contin¬ 
uum level at the emission-line position is estimated through a 
linear fit of the median continuum levels in the continuum re¬ 
gions. This is an important step in the procedure, as the great¬ 
est contribution to uncertainties in the line fluxes stems from 
the determination of the continuum level. Line fluxes are then 

^ At this stage, we only considered as true detections those features 
having simultaneously: (a) a maximum central wavelength shift from 
the expected rest-frame wavelength smaller than 4.5A, and (b) either 
a FWHM within a factor or two, the FWHM of the [On]443727,3729 
doublet, when the doublet was detected, or a FWHM between 2 and 
45A when [O ii]443727,3729 was not detected 


obtained from direct integration of the emission-line profiles 
over the continuum level, which was determined as described 
above. An exception is made for [Nii]46548, Ha, [Nii]46583, 
and [S ii]446717,6731, for which fluxes obtained from Gaussian 
line-fitting were adopted to ensure accurate flux retrieval in cases 
where these sets of lines were blended. 

The preliminary flux measurements were carefully inspected 
for unreliable measurements and/or undetected emission lines 
which are clearly detectable by eye. These checks were done in¬ 
teractively, by visual inspection of the individual spectra and on 
the Gaussian and continuum fits. We also compared our auto¬ 
matic flux measurements with measurements done by hand with 
the SPLOT IRAlQtask for random sets of spectra. The code was 
then modified accordingly until automatic results were in reason¬ 
able agreement (within ~20%) with by-hand flux measurements. 

The uncertainty of the line fluxes was estimated by propagat¬ 
ing the uncertainty of the continuum level emission in the line- 
flux calculation. Emission-line flux measurements with signal- 
to-noise ratios lower than three were rejected and set to zero. 
Typically, relative errors are below ~5% for line fluxes brighter 
than -1.5x10“^^ erg s“^ cm“^ for all emission lines, except for 
[O ii]443727,3729, for which relative errors are larger (-10%) at 
that flux level. Line flux relative errors increase by up to -20% 
at flux levels of - 3 x 10“^^ erg s“^ cm“^ for all lines except for 
[Oii]443727,3729, where errors are typically larger (-25-30%). 

We have empirically checked the quality of our emission¬ 
line measurements by comparing the emission-line ratios, 
[Nii]46548/[Nii]46584 and [0iii]44959/[0iii]45007, to their 
theoretical values. This is done in Fig. The emission-line 
fluxes of both doublets correlate with each other, with correlation 
coefficients of 0.97 and 0.88 for the [Oiii] and [Nii] doublets, 
respectively. The fitted slopes are 0.34 + 0.01 and 0.37 + 0.02 
for [Oiii]44959 vs [Oiii]45007, and [Nii]46548 vs [Nii]46584, 
respectively, which are close to the theoretical value (i.e. 1/3) 
( jOsterbrock & Ferland|2006| ). The dispersion in the [N ii]46548 
vs [N ii]46584 is larger than that in the [O iii] doublet, especially 
for higher values of [N ii]46584/HyS. In general, the galaxies with 
larger deviations from the best linear fits are also those with 
larger error bars. This fact reinforces our flux error estimates, 
which seems to represent the real uncertainties in the measure¬ 
ments well. Data points with the largest deviation from the best 
fit are normally AGNs (especially for the [N ii] doublet). We only 
concentrate on non-AGN galaxies in this paper and, therefore, 
our results will not be affected by these larger uncertainties. 

3.3. Comparison with emission-line fluxes from other 
databases 

We compared our emission-line flux measurements with those 
available in some public databases. The most popular public 
databases with emission-line fluxes are the MPA- JHlfl release 
of spectr um measuremen ts (from SDSS DR7) and the OSS\|^ 
database ( jOh et al.]|2Qll| ). As in our case, these databases also 
used the 7th release of the SDSS spectra. MPA-JHU line fluxes 
of extended sources are normalised to match the photometric 
fibre magnitude in the r-band. We removed this normalisation 
from the MPA-JHU data. 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in 
Astronomy (AURA) under a cooperative agreement with the National 
Science Foundation. 

^ http://www.mpa-garching.mpg.de/SDSS/DR7/ 
http://gem.yonsei.ac.kr/-ksoh/wordpress/ 
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Fig. 1. Comparison of the measured emission-line fluxes of the two lines in doublets [N ii]/l/l6548,6584 {left) and [O iii]/l/l4959,5007 {right) 
normalised to the emission-line flux. Purple points indicate galaxies classifled as AGN (Sect.j^. The dashed green line shows the best linear 
fit to the data, which yields a slope of 0.37 + 0.02 and 0.34 + 0.01 for [N ii]/l6548/Hf5 vs [N ii]/l6584/Hf5, and [O iii]/l4959/Hf5 vs [O iii]/l5007/Hf5, 
respectively (see text in Sect. |3.2| for details). 
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Fig. 2. Ratio of emission-line fluxes available in public databases to the fluxes measured in this work (FIux^vk), for the OSSY (green circles) 
and MPA-JFIU (pink squares) databases, as a function of the emission-line fluxes measured in this work for the [O ii]/l/l3726,3729 (top left), 
[O III ] /t50 07 (top right). Her (bottom left), and [Nii]/16583 (bottom right) emission-lines. All fluxes are corrected for Galactic extinction (see 
Sect. |3.3| for details). 


Figure shows the comparison of the [Oii]/l/13727,3729, 
[Oiii]T5007, Her, and [Nii]/16583 emission-line fluxes, mea¬ 
sured by these databases, with our measurements from the spec¬ 
tra of our galaxy sample (Sect. [^. The three sets of data rep¬ 
resented in Fig. have been corrected for foreground Galac¬ 
tic reddening (but not for internal extinction). Our emission-line 
fluxes match fairly closely those from the OSSY and MPA-JHU 
databases. In general the fluxes match within ~20% for lines 
brighter than erg s“^ cm“^, except for [O ii]/l/13727,3729, 

where differences are within ~30%. For weaker lines, the scatter 
increases considerably, as a result of larger relative errors in the 
measurements. 

For emission lines brighter than >3x10“^^ erg s“^ cm“^, our 
fluxes match with the OSSY and MPA-JHU measurements to 
within ~4-6%, except for [Nii]/16583 and Ha as measured by 
OSSY, which are, on average, ~8% higher than our measure¬ 
ments. Our measurements of [Nii]/I6583 and Ha are also sys¬ 
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tematically lower than those from MPA-JHU, but only by ~4% 
on average. For [Oii]/l/l3727,3729 and [O iii]/l5007, however, 
our measurements are systematically ~5% higher than MPA- 
JHU measurements, and are a better match with the OSSY mea¬ 
surements (within ~2-3%). These systematic differences can, 
in part (but not completely), be explained by the use of differ¬ 
ent extinction laws for the correction of the foreground Galac¬ 
tic extinction. Our data was corrected using the |Schlegel et af 


( 1998) maps and the|Cardelli et al.| ( |1989[ ) extinction curve (see 
|Cid Fernandes et al.|20()5|). The MPA-JHU data have been cor- 
rected using |0’Donnell| ( [1994] ). Differences of approximately 
~3% at [Oii]/l/13727,3729 are expected for typical values of 
Ay~ 1-1.5 mag, given the differences in the extinction curves. 
The extinction curve used by OSSY authors is not specified. 


In any case, differences between our observed fluxes and the 
ones available in the OSSY and MPA-JHU databases are small 
and comparable to measurement errors, even taking into account 
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that the stellar emission subtraction, the flux measurement pro¬ 
cedures, and the foreground Galactic reddening correction differ 
in the three sets of data. 


4. Nuclear classification and removai of AGN 

Our initial galaxy sample (see Sect. ^ includes active and nor¬ 
mal galaxies. We used standard diagnostic diagrams based on 
bright emission-line ratios (also known as BPT diagrams, after 
[Baldwin et al.||198T] ) to separate the galaxies where the dom¬ 
inant mechanism of ionisation is UV radiation from massive 
young stars (i.e. star-forming galaxies) from those galaxies being 
ionised by an active galactic nucleus (AGN). In Figs. and [Tb] 
we show the diagnostic diagrams based on our measured (non¬ 
extinction correctec0) ratios [Oiii]/15007/HyS vs [Nii]/I6584/H(r, 
for the barred and unbarred sub-samples of galaxies. 

The solid line separates star-forming galaxies from AGN ac- 


cording to the predictions of the photoionisation models of Kew- 
jley et al. (|2001|), while the dash ed line sets this separation empir- 


ically ( jKauffmann et al.pOOS^. Both lines do not match , and the 
number of AGN predicted by jKauffmann et al.| ( |20Q3a| ) is larger 
than the theoretical prediction. We refer to the galaxies located 
betwe en the two lines (i.e . classified as ’’star-formi ng" accord¬ 
ing tojKewley et al.|2001| and as AGN according to jKauffmann 


jet al.|2003a| ) as transition objects. These are also known as com¬ 
posite objects and might be a population of objects where ionisa¬ 
tion is partly produced by recent star-formation and partly due to 
an AGN ( jKewley et al. [20061 ), but they might al so include pho- 
toionised nuclei with high nitrogen abundance (jPerez-Montero 


& ContinT||2009| ). For this reason we used the jKewley et al. 


^ 20Qlj ) criteria to remove AGN from our sample. The horizon 


tal dot-dashed line separates the Seyfert region from the Low- 
Ionisation Nuclear-line Regions area (LINERs; jKewley et al. 
[2006] ). 

According to this classification into star-forming, transition 
and AGN galaxies, the percentage of AGN in the sub-samples 
of barred and unbarred galaxies is 31% and 18%, respectively. 
In other words, there are ~ 1.7 times more AGN in barred 
than in unbarred galaxies. This factor remains when we consider 
the jKauffmann et ( j2Q03aj ) criteria, i.e. AGN (according to 
[Kewley et al.||200l] rplus transition objects are 1.7 times more 
frequent in the sub-sam^e of barred galaxies than in the un¬ 
barred galaxy sub-sampl^ The fraction of AGN in barred galax¬ 
ies increases when we consider galaxies with low-mass bulges. 
AGN are about twice as frequent in barred galaxies as in un¬ 
barred galaxies, if we only consider bulges with mass lower than 
lO^^-^Mo. 

Below, we consider all the galaxies of the sample, except 
those classified as AGN. We will refer to them as non-AGN 
galaxies, and this sample includes star-forming galaxies accord¬ 
ing to the criteria of |Kauffmann et al.| ( j2Q03aj ), transition objects, 
and unclassified galaxies (normally due to the non-detection of 


^ The [O III] 25007/1^ vs [N ii]26584/Hq' is the least reddening de¬ 
pendent BPT plot. After internal extinction correction as described in 
Sect. 1^ we have checked that our classification was correct. BPT dia¬ 
grams shown in Sect.[^are based on extinction-corrected emission-line 
fl uxes. _ 

^ Coelho & Gadotti (|201l|) u sed t he AGN cata¬ 
logue from jKauffmann ^ ah] j2003aj http://www.mpa- 
garching.mpg.de/SDSS/DR4/Data/agncatalogue.html to separate 
AGN galaxies from the whole sample. Their figures differ from 
the ones used in this work, in which classification comes from our 
emission-line measurements. However, in both cases AGN are ~ 1.7 
times more frequent in barred than in unbarred galaxies. 


Table 1. Classification of the galaxy sample. 


Type 

whole 

sample 

barred 

subsample 

unbarred 

subsample 

Star-forming 

139 

43 

96 

Transition 

113 

66 

47 

LINER 

109 

61 

48 

Seyfert 

28 

17 

11 

Unclassified 

186 

64 

122 


Notes. Number of active and non-active galaxies in the global sample 
and in the sub-samples of barred and unbarred galaxies, according to 
the BTP diagrams in Figs.and [^ (see Sect, j^ for details). 



Fig. 4. Histograms showing the distribution of decimal logarithm of the 
total galaxy stellar mass, redshift, logarithm of bulge mass and loga¬ 
rithm of the B/D light ratio in the i band for barred (hatched yellow) 
and unbarred (purple) galaxies separately, for all the non-active galax¬ 
ies of the sample. The number of galaxies in each sub-sample (N^ and 
Nf/) is indicated, together with the two-sample Anderson-Darling test 
P-values (AD, expressed in %). The yellow and purple circles in the 
upper side of the panels indicate the median value of the barred and 
unbarred distributions, respectively. The horizontal error bar covers the 
95% confidence interval (estimated as 1.57 x IQR/ ^^N, where IQR is 
the interquartile range, or 1st quartile subtracted from the 3rd quartile, 
and N the number of data points) for the corresponding median value. 


relevant emission lines involved in the BPT diagram). Our effec¬ 
tive sample of non-AGN galaxies has 173 and 265 barred and 
unbarred galaxies, respectively. 


4.1. The characterisation of non-AGN and star-forming 
sub-samples 

Figure shows a comparison between the total galaxy stellar 
mass, redshift, bulge mass, and bulge-to-disc (B/D) light ratio 
(/-band) distributions for the sub-samples of barred and unbarred 
galaxies. The total stellar masses are taken from |Kauffmann et aL 
( 2003bj ) and are based on fitting to stellar spectral features (the 
4000A break and the Baimer absorption line index H^a). The 
bulge masses come from jCoelho & Gadottij ( j2011| ) (based on 
bulge luminosity and mass-to-light ratio in the / band) and the 
B/D light ratio comes from morphological decomposition per¬ 
formed by [Gadottij ( [2009 j ). The total galaxy stellar mass and red¬ 
shift distributions are similar for barred and unbarred galaxies. 
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Fig. 3. Diagnostic diagrams to classify emission-line galaxies according to their dominant ionisation mechanism for the barred (leff) and unbarred 
(right) galaxies of the sample. The solid and dash lines separate the region of star-forming galaxies from the region of AGN according to lKewle}^ 
|et al.| ( [200T| and|Kauffmann et al.| ( [20()3a| ), respectively. The horizontal dot-dashed line separates the AGN region into Seyfert and LINER galaxy 
regions ( [Kewley et al.|2006| ). 



loglM./Mg) Redshift 

1.6 

1.4 

1.2 

^ 1.0 
0) 

g- 0.8 

OJ 

it 0.6 
0.4 
0.2 

Fig. 5. Same as Fig.[^but only for gala xies classified as "star-forming" 
according to |Kauffmann et al!] ( |2003a| ). N^ and Nj/ are the number of 
barred and unbarred galaxies, respectively, and AD is the two-sample 
Anderson-Darling test P-value (expressed in %). 



To confirm this statistically we used the ^-sample Anderson- 

Darling (hereinaf ter A-D) test (see |Scholz & Stephens||1987 

■ . 


Stephens||1974||^| It tests the null hypothesis that the two sam¬ 
ples are drawn from the same parent distribution. The output of 


^ The ^-sample Ander son-Darling test is based on the Anderson- 
Darling rank statistics ( [Anderson & Darling||1954] | for testing ho¬ 
mogeneity of samples with possibly different sample sizes and un¬ 
specified distributions. The Anderson-Darling test is, in turn, a mod- 
ifieation of the Kolmogorov-Smimov test. The /:-sample version of 
the test is more recommended than ^-sample Kolmogorov-Smirnov 
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the A-D test is a P-value or significance level at which the null 
hypothesis can be rejected. For the total stellar mass and red- 
shift, we find P-values of ~37% and 54%, respectively. These 
values are much greater than the 5% threshold value normally 
adopted below which the result is statistically significant and the 
null hypothesis can be rejected. The distributions of bulge mass 
and B/D light ratio are, however, different between barred and 
unbarred galaxies when all non-AGN galaxies are considered. 
Median values and distribution shapes are different and the A-D 
tests confirms this difference, with P-values well below 0.1%. 
The bulges of barred galaxies are less massive than in unbarred 
galaxies and B/D light ratios are also lower for barred galaxies, 
even though both sub-samples of barred and unbarred galaxies 
have the same total stellar mass distributions. This has already 
been pointed out by jCoelho & Gadotti|p011| ). 

If we consider the sub-sample containing only galaxies clas¬ 
sified as star-forming, Fig.j^and the A-D test results indicate that 
the distributions of total stellar mass, redshift, bulge mass, and 
B/D light ratio are no different for barred and unbarred galaxies, 
since P-values are higher than 5% in ah cases. 


5. Internal extinction 

The emission-line fluxes were corrected for internal extinction 
from the Baimer decrement using the Htr to H/5 emission-line 
flux ratio. In the absence of internal dust extinction, the Ha to HjS 
line flux ratio, (pHa/^Hfi)intr is equal to 2.86, for typical electron 
temperatures and densities in star-forming regions (T^ -10000 K 
and xie -100 cm“^), in the case B of the recombination theory 

. The presence of dust in the inter¬ 
line ratio as a result of differential 
extinction. The amount of attenuation can be quantified with the 

for small sample sizes (e.g. |Hou et aL]|2009j ). See also discussion at 
https://asaip.psu.edu/Articles/beware-the-kolmogorov-smirnov-test. 


( Qsterbrock & Ferland|2006| ) 
stellar medium increases this 
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Fig. 6. Left. Histogram of the values of the B aimer extinction at the Yip 
emission line, ciYip), for all galaxies of the sample with ciYlp) >0.0, 
except those classified as AGN as explained in Sect. Right. Same 
as for the left, except for barred (hatched yellow) and unbarred (purple) 
galaxies separately. The number of galaxies in each sub-sample, and the 
P-value from the A-D test for the two distributions, are shown below the 
figure legend. 


parameter c(HyS), the internal extinction at the H/5 emission line 

{FhhUs = ( 1 ) 


where the sub-indexes "obs" and "intr" stand for observed and 
intrinsic fluxes, respectively. The Yla to Yip emission-line ratio 
together with an assumed extinction curve /(T) can be used to 
obtain c(HyS) 


cim = 


-1 

/(Ha) 


log 






mimv 


( 2 ) 


where /(Ha) is the reddening function at Ha normalised to Yip 
(i.e. f( HP) = 0). We emp loyed the|Seato^([1979|) reddening law 


with the 


Howarth 


emp loyed the |Seato] ^( [1979| ) 
(|1983) parametrisatioijj assun 


assuming Ry = 3.1. 


With this combination of reddening law and Ry, the extinction 
in magnitudes at Ha is Aua = 1.515 xc(HyS). The intrinsic Ha to 
Yip line flux ratio was assumed to be 2.86. We are aware that in 
the central regions of galaxies the metallicity might be oversolar, 
with Tg lower than 10000 K, and, therefore, this ratio may be 
higher than 2.9. Assuming an increased value of the Ha to Yip 
line flux ratio would decrease all c(Hy5) values by 0.02, 0.06, and 
0.1 for (^Hal^Hi 5 )intr, cqual to 2.9, 3.0, and 3.1, respectively. 

The histogram with the values of the B aimer extinction 
c(HyS) for all galaxies of the sample (except AGN) is in the left- 
hand panel of Fig.[^ Values occur in the range 0-1.4, being the 
median value of the distribution (0.47+0.03). We assign an inter¬ 
nal extinction of zero when the observed Baimer ratio Vnal^Hp 
is lower than the theoretical value. These objects are not consid¬ 
ered in Fig.[^ The distribution obtained is very similar in shape 
and values to that found by [Stasinska et al.| ( |2004| ) for 10854 
spectra from the First Data Release of the SDSS. 

Dust affects both the stellar continuum emission and the neb¬ 
ular emission. Figure [7] shows the V-band extinction in magni¬ 
tudes derived from the gas emission lines (Ayg) versus the stellar 
continuum V-band extinction (Ay^^), as derived from the spec¬ 
tral synthesis fitt ing with STARLIGH T (see Sect.j^. We should 
point out that the |Cardelli et al.| ( [l989| ) extinction curve was used 


The relation of /(T) with the function X(x) of the 


parametrisation is /(T) = “ F where x 


Howarth 




l/A (A in microns). 


with xh /3 the corresponding value for the spectral line. In this 
parametrisation, the extinction in magnitudes at a given wavelength A 
is A^ = X(x) E(B - V). 



Fig. 7. V-band extinction in magnitudes, Ayg, as obtained from the 
Her to Yip B aimer decrement, versus the stellar continuum V-band 
extinction (AyJ, as derived from the spectral synthesis-fitting with 
STARLIGHT for all non-AGN galaxies. The straight sol id and dashed 
lines s how the relations betwee n Ayg and Ays found by |Calzetti et al.| 
( |2000| ) and |Kreckel et al.|p013| ), respectively. 


in STARLIGHT, and we used |Seaton| \1919) for the gas com¬ 
ponent. However, both extinction curves give roughly the same 
extinction in the V-band (Seaton gives on average a higher ex¬ 
tinction rate by -0.8% in that spectral range). The dispersion is 
high, but it is clear from the figure that both quantities are cor¬ 
related, and that the extinction derived from emission lines is 
typically about twice the extinction derived from the observed 
stellar continuum. This is in agreement with other author results, 
such as |Calzetti et al.| ( [2000| ) and |Kreckel et al.| ( |2013| ), who find 


Ay^=(0.44 ± 0.03) Ay^ and Ay^=(0.470 ± 0.006) Ay^ for Hii 
regions and starburst galaxies, respectively. These relations are 
overplotted in Fig. [7] and are compatible with our data in the 
range of extinction observed (Ay^ from 0 to 3 mag). 


All measured emission lines were corrected for internal dust 
extinction using the derived values for c(HyS). The corrected 
emission-line fluxes have been compared with the correspond¬ 
ing values provided in one of the most popular public databases: 
OSSY dOh et al.|[20TT] ). Figure [8] shows the comparison of 
extinction-corrected fluxes for some of the brightest emission 
lines. As can be seen, OSSY fluxes are typically larger than our 
extinction-corrected line fluxes by a factor of about 1.6 (-1.8 
for [On]). Also, the OSSY colour excess tabulated values for 
the gas and the stellar continuum (called E(B-V)_gas and E(B- 
V)_star, respectively) do not maintain the relation mentioned 
above (Fig. [7]). However, we would like to point out here that 
our observed (uncorrected) emission-line fluxes correspond ex¬ 
tensively with OSSY quoted observed fluxes (see Sect. |3. 3 [ and 
Fig-Hf. 

A table containing extinction-corrected line fluxes measured 
in this paper is available at the CDS. The table contains the 
following information: Column 1 gives the SDSS plate used 
to collect the spectrum; Column 2 the modified Julian Date of 
the observation night (mjd); Column 3 the SDSS fibre number; 
Columns 4 to 25 the extinction-corrected emission-line fluxes 
normalised to H/5=100 and associated uncertainties; Columns 
25 and 27 give the Yip flux in units of 10“^^erg s“^ cm“^ and 
the Baimer extinction at Yip, respectively; and Columns 26 and 
28 give their corresponding uncertainties. 
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Fig. 8. Ratio of internal extinction corrected emission-line fluxes from the OSSY database to those measured in this work, as a funetion of the 
latter for the [Oii]/l/13726,3729 (top left), [Oiii]/I5007 (top right). Her (bottom left) and [N ii]/l6583 (bottom right) emission lines. All fluxes have 
been corrected for Galactic extinction and internal extinction (see Sect.[^for details). 


5.1. The distribution of c(Hj3) for barred and unbarred 
galaxies 

The right-hand panel of Fig. shows the distribution of c(U/3) 
for barred and unbarred galaxies separately. The median value 
of c(Uj3) is only marginally larger in barred galaxies (0.50+0.04) 
than in unbarred (0.44+0.03), but the A-D test indicates that the 
two distributions differ, with a P-value of 0.8%, which is well 
below 5%. We find the same result when we only consider star¬ 
forming galaxies (see Table [^. 

Previous work on the Baimer extinctio n at Hj3 on large SDSS 
data samples (e.g. [Stasinska et al.|[20Q4] ) shows that c(HyS) de¬ 
creases from early- to late-type galaxies. Additionally, c(Hj3) 
also depends on total galaxy luminosity, being higher for the 
most luminous galaxies (and probably more massive). Our sub¬ 
sample of barred non-AGN galaxies is similar to the unbarred 
sub-sample in terms of total stellar mass (Fig.|^, but it is biased 
towards lower bulge mass and lower B/D light ratios compared 
with the sub-sample of unbarred galaxies. If we assume that the 
results obtained by [Stasinska et al.| ( |2004| ) for integrated spectral 
properties of galaxies are applicable to the inner kiloparsec of 
galaxies, we would predict a lower average of c(Uj3) in barred 
than in unbarred galaxies. However we find just the opposite, al¬ 
though differences in median values are only marginal (within 
error bars). 

The sub-samples of barred and unbarred star-forming galax¬ 
ies (Fig.[^ are equivalent in B/D light ratios, and bulge mass and 
their distributions of c(HyS) are also different (P- value 0.6%), be¬ 
ing larger in barred than in unbarred galaxies by ~ 0.07 dex, ap¬ 
proximately the same difference as between the non-AGN barred 
and unbarred sub-samples. A larger central dust extinction might 
be related to a larger dust mass surface density ( jKreckel et al.j 
|2013| ) and, therefore, the observed difference in c(Hj3), although 
marginal, could be due to different central dust mass concentra¬ 
tions possibly resulting from the transfer of material towards the 
galaxy centres, induced by the bars. 
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6. Central oxygen abundance and N/0 abundance 
ratio 

The faint auroral emission lines are not generally detected in our 
optical spectra and, consequently, the electron temperature (Te) 
cannot be calculated. Therefore, it is not possible to determine 
abundances from the -based or direct method and the oxygen 
abundance can only be inferred through empirical or theoreti¬ 
cal calibratio ns of nebular strong-line flux ratios, the strong-line 
methods (see jPerez-Montero & Diaz|2005[|L6pez-Sanchez et al. 


|2012| for in-depth discussions on different methods). We have 
used two of the most widely used tracers of gas oxygen abun¬ 
dance, R 23 = ([Oiii]TT4959,5007+[Oii]TT3727)/HyS, and N2 
= log([Nii]/16583/HQf). The R 23 parameter was first proposed 
by jPagel et al.| ( |1979] ), but there are many diff erent empirical and 
theoretical calibrations of this parameter (e.g.|Pilyugin|200f||Pi-| 


Kobulnicky & 


lyugin & Thuan 2005 [|McGaugh|l991] [Kewley & Dopita|2002t 


■Cewleyl 


_ ||2004|). The N2 was first proposed by 

Storchi-Bergmann et al.j ( j 199^ and has been extensively used, 
especially for intermediate and high redshift studies. This is be¬ 
cause the lines involved are readily observable using the current 
generation of telescopes, and the line ratio is insensitive to red¬ 
dening corrections due to the small wavelength separation be¬ 
tween the two involved lines. Some of the most common cali- 
brations of N2 include jPenicolo et al.j (|2Q02j), [Pettini & Pagelj 


( 2004| ), [Nagao et al.] ( |2006| ) or jPerez-Montero & Contini| ( |2009| ). 

Our aim is to investigate potential differences between barred 
and unbarred galaxies, and therefore what follows only com¬ 
pares strong-line oxygen abundance indicators, without apply¬ 
ing any of the existing calibrations. Figure shows the com¬ 
parison of the distribution of values for the R 23 and N2 param¬ 
eters for barred and unbarred non-AGN galaxies. Barred and 
unbarred galaxies seem to be indistinguishable in terms of the 
median value and distribution of the R 23 parameter. There is, 
however, a striking difference in N2 where barred galaxies have 
enhanced log([Nii]/16583/HQ') compared with unbarred galax- 
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Table 2. Statistics for barred and unbarred non-AGN and star-forming galaxies. 




Non-AGN 




Star-forming 



P-value'^ 

barred 

median 

N 

unbarred 

median 

N 

P-value^ 

barred 

median 

N 

unbarred 

median 

N 

dm 

0.8% 

0.50+0.04 

104 

0.44+0.03 

140 

0.6% 

0.53+0.04 

43 

0.46+0.04 

96 

log R 23 

64.0% 

0.31+0.04 

105 

0.31+0.04 

130 

0.3% 

0.19+0.04 

41 

0.22+0.04 

86 

N2 

0.0009% 

-0.28+0.02 

no 

-0.38+0.02 

152 

0.02% 

-0.37+0.03 

43 

-0.43+0.01 

96 

log '^SFR 

0.07% 

-1.1+0.1 

112 

-1.3+0.1 

155 

0.1% 

-0.9+0.1 

43 

-1.2+0.1 

96 

[Sii]i6717/[Sii]^6731 

0.05% 

1.25+0.03 

104 

1.31+0.03 

127 

0.7% 

1.28+0.03 

40 

1.33+0.03 

84 

log (N/O) 

0.01% 

-0.49+0.03 

105 

-0.58+0.03 

143 

0.03% 

-0.49+0.04 

40 

-0.58+0.03 

91 

12+log (0/H) 

4.3% 

8.68+0.01 

107 

8.69+0.01 

151 

7.4% 

8.69+0.01 

40 

8.70+0.01 

92 

logU 

5.3% 

-3.08+0.02 

107 

-3.14+0.02 

151 

2.5% 

-3.09+0.03 

40 

-3.18+0.02 

92 


Notes. Median values and errors (95% confidence interval of the median) of the distributions of: c{]Ap), R 23 and N2 gas-phase metallicity indicators, 
star formation rate per unit area, [S ii]/l6717/[S ii]/16731 line ratio, oxygen abundance and N/0 ratio, and logarithm of the ionisation parameter. 
All figures refer to the centres of barred and unbarred galaxies in the sub-samples of non-AGN and pure star-forming galaxies (Sect.|^. N is the 
number of galaxies in each sub-sample. The P-values of the ^-sample A-D test for the comparison of distributions for barred and unbarred galaxies 
are shown for all cases. 

(«) A -D test P-value or approximate significance level at which the null hypothesis (the two samples are drawn from the same population) can be 
rejected. Usually, significance levels lower than 5% are requested to reject the null hypothesis. 



empirical calibration by |Perez-Montero & ContinT (2009), which 
includes a correction to account for the N/0 ratio. We have esti¬ 
mated log (N/0) from the empirical calibration of the N2S2 ra¬ 
tio (=log([N ii]/f6583/[S ii]AA6111, 6731) from the same authors. 
Afterwards, 12 -h log(0/H) is calculated. We have also esti¬ 
mated log (N/0), 12 -hlog(0/H) and the ionisatio n parameteip] U, 


from the new grids of photoionisation models by |Perez-Montero 
( |2014| ). The results are shown in Table and in Fig. |10[ The 


-0.2 0.0 0.2 0.4 0.6 

log R23 


Fig. 9. Histograms of the values of the oxygen abundance indicators 
P 23 =([Oiii]+[Oii])/Hf 5 (left) and N2 = log([Nii]T6583/HQf) (right) for 
all non-AGN barred and unbarred galaxies in the sample. The number 
of objects in each sub-sample and the P-values from the A-D test are 
shown in each panel. 


ies by 0.10 dex (see Table |^, with both displaying significantly 
different distributions according to the A-D test (P-value be¬ 
low 0.001%). The two distributions of N2 are also different if 
we only consider star-forming galaxies (P-value 0.002%), and 
barred galaxies also have a larger average N2 than unbarred 
galaxies. However, the barred and unbarred difference in median 
values shortens to -0.06 dex in the star-forming sub-sample, 
which is, in any case, much larger than the uncertainties in the 
median values. 

It is well known that there is a positive mass-metallicity rela¬ 
tion for spirals (e.g. |Tremonti et al.|2004[|Bothwell et al.|2013| ). 
We will analyse in detail the dependence of N2 and R 23 with to¬ 
tal stellar mass in Sect. [T^ However, it should be noted that the 
mass distributions of our barred and unbarred sub-samples are 
identical and therefore the observed difference in N2 might indi¬ 
cate a real and strong physical difference in the centres of barred 
galaxies compared with unbarred galaxies. 

Given that there are no differences between barred and un¬ 
barred galaxies in terms of the R 23 parameter, any empirical cal¬ 
ibration of this parameter yields no difference in oxygen abun¬ 
dance between barred and unbarred galaxies. The opposite hap¬ 
pens with calibrations that make use of N2, except for the N2 


barred and unbarred sub-samples of galaxies have the same aver¬ 
age oxygen abundance within errors, 12-rlog(0/H) = 8.69+0.01, 
and the distributions are indistinguishable from each other: we 
obtain A-D P-values of 4.3% and 7.4% when we consider non- 
AGN or star-form ing galaxies, respecti vely. This agrees with re¬ 
cent results from |Cacho et al.| ( |2014| ) for a similar sample of 


barred and unbarred nearby disc galaxies from the Nair & Abra- 
|ham (|2010|) catalogue. H owever, our result apparently conflicts 
with Ellison et~ar] ( |2011| ), who report a larger oxygen abundance 
(by -0.06 dex) in barred galaxies compared to unbarred galax- 
ies, in relation t o a galaxy sample also extracted from |Nair &| 
[Abraham] ( |2010| ). In Sect. [T^ we discuss these results and try to 
explain the possible sources of discrepancy between their results 
and ours. 

However, as expected from the N2 values, we find a statis¬ 
tically significant and interesting difference in the central N/0 
abundance ratio found in barred and unbarred galaxies. The me¬ 
dian log (N/0) is 0.09 dex larger in barred than in unbarred 
galaxies, for both the non-AGN and the star-forming samples. 
This difference is three times larger than the uncertainties in the 
median values (-0.03 dex). The A-D test confirms that the barred 
and unbarred distributions are different (mid-panel in Fig. [T^, 
with P-values of 0.01% and 0.03% for the non-AGN and the 
star-forming samples, respectively. To our knowledge, this is the 


^ Regarding the ionisation parameter values derived from the grid of 
photoionisation models from jPerez-Montero] ( |2014| ), it is important to 
stress that, in absence of emission-line ratios sensitive to the electron 
temperature, an empirical law between log U and 12+log(0/H) is as¬ 
sumed (i.e. lower U for higher 0/H). Nevertheless, in the grid a certain 
range of variation in log U is allowed for each 0/H value that exceeds 
the expected and reported log U variations in our sample of studied ob¬ 
jects. As such, we think that the resulting log U values give an accurate 
idea of the variations. 
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12 + log {0 /H)pmi4 




Fig. 10. Histograms of the oxygen abundance (top), N/0 ratio (mid- 
dle) and ionisati on parameter (U, bottom), as calculated with the |Perez-] 
|Montero| ( |2014| ) method for all non-AGN barred and unbarred galaxies 
in the sample. The number of objects in each sub-sample and the P- 
values from the A-D test are shown below the legend. Only non-AGN 
galaxies are considered. 


most obvious and largest difference in physical properties of the 
gas so far observed between barred and unbarred galaxies. 


7. Star formation rate 

We have estimated the star formation rate (SFR) in the galaxy 
centres from the Htr extinction-correct ed emission within the 3" 
SDSS fibres and the Kennicutt ( |l998 ) conversion factoipl The 
distances to the galaxies have been estimated from the redshift 
given by SDSS and a Hubble constant of 72 km s“^ Mpc“^. 


The stellar component had already been subtracted from the origi¬ 
nal SDSS spectra and, therefore, we do not need to correct for stellar 
absorption under the Baimer lines at this point. 
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Fig. 11. Comparative histograms of star formation rate per unit area for 
all non-AGN barred (yellow hatched) and unbarred (purple) galaxies of 
the sample. 


Given the redshift range of the galaxy sample, the SDSS 3" 
diameter fibre corresponds to a projected size ranging from -1.2 
to 4.2 kpc. Although the redshift range is small and the redshift 
distributions of barred and unbarred galaxies are similar, we have 
derived star formation rates per unit area, ^sfr, to average out 
possible dependencies on galaxy distance. 

The distributions of the logarithm pR for barred and un¬ 
barred galaxies are shown in Fig. [TT] Both distributions are dif¬ 
ferent according to the A-D test (R-value < 0.1%), with barred 
galaxies showing marginally larger star formation rate per unit 
area. This difference is slightly larger when considering only 
star-forming galaxies, where barred galaxies show a median 
log ^SFR which is -0.3 dex larger than unbarred galaxies. A 
lar ger SFR in the cent res of barred galaxies was also reported 
by [Ellison et al. |( |2Qll| ) for galaxies of similar total stellar mass 


in a similar study, and by other authors (see e.g. Hummel et al. 


Ho et al.||1997bl|Wang et al.||2012t|Oh et al.||20121|Zhou 

joBF 


1990 


et al. 


8. Electron density 

The [Sii]/16717/[Sii]/16731 line ratio is sensitive to the electron 
density (Ng). Nearly 8% of all galaxies in our sample have a 
[S ii]/16717/[S ii]/16731 line ratio which is above the low-density 
theoretical limit (=1-43 , implying -lcm“^, for T^^IOOOOK; 
Shaw & Dufour|p~994| ). Most of these targets have the largest 
[S ii]/16717/[S ii]/16731 relative errors (-20-40%). Here, we only 
consider galaxies where their [S ii]/16717/[S ii]/16731 line ratio 
value is still compatible with being below the low-density the¬ 
oretical limit, taking into account their corresponding error bars 
( jLopez-Hemandez et al.|2013] ). The median value of this line ra¬ 
tio in the centres of the non-AGN galaxy sample is 1.29+0.02 
(1.30+0.02 for pure SF galaxies), which corresponds to N^ - 
100 - 150 cm“^. This is in agreement with typical values found 
in centres of galaxies measured previously by [Kennicutt et al. 
( 1989| ) and Ho et al.[ ( [1997a| ). In Fig. [T^ we compare the dis¬ 
tributions [Sii]/l6717/[Sii]/16731 for barred and unbarred non- 
AGN galaxies. The median value is marginally lower for barred 
(1.25+0.03) than for unbarred galaxies (1.31+0.03), indicating 
that, on average, barred galaxies tend to posses a larger central 
electron density than unbarred galaxies. The A-D test yields a 
very low R-value (0.05%), which indicates that both distribu- 
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Fig. 12. Comparative histograms of values of the [S ii]/l6717/[S ii]/l6731 
line ratio for barred (yellow hatched) and unbarred (purple) galaxies 
separately. Only non-AGN galaxies are considered. 


tions are different. Distributions are also different in the sub¬ 
sample of pure star-forming galaxies (P-value of 0.7%, see Ta- 
ble[|). 

9. Dependence of central differences on Hubble 
type 

The amount of gas driven by a bar towards the galaxy centre 
might depend upon, amongst others, the availability of gas in 
galaxy disc and on the gravitational torque created by the bar¬ 
like mass distribution. Gas content varies across the Hubble se¬ 
quence, but also bars in early-type galaxies tend to be longer and 
stronger than in late-type spirals ( Erwin|[2005t [Chapelon et alT 


1999). Our sub-samples of non-AGN galaxies are not similar in 


terms of B/D flux ratio and bulge mass (Fig. |^, which might 
indicate that both sub-samples are different in the relative popu¬ 
lation of morphological types. We next explore whether any of 
the differences between barred and unbarred galaxies found in 
previous sections could be due to a bias towards earlier types in 
the sub-sample of barred galaxies. 

The fact that most of the differences between barred and un¬ 
barred galaxies of the non-AGN sub-sample are also observed 
in the pure star-forming sub-sample, where the distribution of 
bulge mass and bulge-to-disc flux ratio are equivalent for barred 
and unbarred galaxies (Fig. [^, seems to indicate that, even if a 
bias were present, bars tend to produce higher dust concentra¬ 
tions, larger N/0 ratios, star formation, and electron densities in 
the centre of galaxies. A deeper analysis is, in any case, neces¬ 
sary to clarify the nature of these differences. 

Morphological classification is only available for -46% of 
the total galaxy sample, i.e. 245 galaxies, of which 113 are 
barred and 132 unbarred. This classification comes from the lNairl 
& Abraham] ( |2010| ) catalogue. However, morphological decom¬ 


position of the galaxy sam ple has already been performed using 
the code BUDDA v2.l[^ ^de Souza et al.|[2004 Gadotti|[2008 


|2009| ) that is able to fit up to four different galactic components 
(a bulge with a Sersic profile, a single or double exponential disc, 
a Sersic bar, and a Moffat central component) to a galaxy image. 
The morphological decomposition of the galaxies from t his sam- 
ple was performed over SDSS g, r, and /-band images ( |Gadotti| 

^^ http://www.sc.eso.org/-dgadotti/budda.html 


2009| ) and the B/D flux ratios for the three photometric bands 
were therefore obtained. 

The behaviour of the B/D flux ratio along the Hubble se¬ 
quence is well known. However, this relation has a high dis¬ 
persion and varies significantly, depending on the function used 
to fit the bulge vs Sersic on the photometric band 
or on the galaxy inclination and dust-extinction (see [Graham & 


I Worley [2008 1 and references therein). We have analysed the de¬ 
pendence of the logarithm of the g, r, and /-band B/D light ra¬ 
tios obtained with BUDDA ( |Gadotti|[2009 ) with the T-type by 
[Nair & Abrahamj ( |2010| ) for the galaxies in common between 
both samples. From these relation s and additional inf ormation 
from the Galaxy Zoo catalogue ( j Willett et al. 2013| , we have 
separated ah galaxies in our sample into early- (T-type<2) and 
late-types (T-type>2). For more information on this separation, 
please see Appendix |A| 

The bulge stellar light and, therefore, its mass has been 
claimed as a key parameter that might distinguish galaxies across 
the Hubble se quence (e.g. [Meisels & Q striker] 1984 [ [Graham ^ 
[Worley[2008 1 and references therein). In addition to the B/D light 
ratios we have also used the bulge mass ( [Coelho & Gadotti[201 1[ ) 
to explore whether the observed differences between barred and 
unbarred galaxies are dependent on Hubble types. We have sep¬ 
arated ah non-AGN galaxies into two groups according to their 
bulge mass. The dividing bulge mass (10^-^ M©) has been se¬ 
lected to ensure that the total stellar galaxy mass distributions 
for barred and unbarred galaxies are equivalent in both bulge 
mass bins (with A-D F-values greater than 10%). 

Figure and Table show a comparison of central proper¬ 
ties between barred and unbarred galaxies separately for early- 
and late-type disc galaxies (using B/D light ratios, as explained 
above). The same comparison is made in Fig. and Table 
for galaxies with bulge mass below lO^-^M© and bulge mass 
larger than or equal to 10^-^ M©. It can be seen that for galax¬ 
ies with later types (or less massive bulges), the distributions 
of ah parameters except log R23 are different between barred 
and unbarred galaxies, with A-D F-values lower than 4% in ah 
cases. For earlier type galaxies results depend on whether we 
split galaxies attending to B/D light ratios or bulge mass. Barred 
and unbarred earlier-type galaxies (according to B/D) only show 
different distributions for N2 and the [S ii]46717/[S ii]46731 
emission-line ratio, while for galaxies with bulge mass above 
lO^-^M©, in addition to those parameters, we also hnd different 
distributions for IIsfr and log (N/0). 

Barred galaxies tend to have, on average, larger central 
B aimer H/3 extinction, larger values of the N2 metahicity indi¬ 
cator, SFR per unit area, electron densities and N/0 abundance 
ratio than unbarred galaxies in ah sub-samples. However, taking 
into account the uncertainties in the median values, we cannot 
claim real differences in average values between these param¬ 
eters in barred and unbarred galaxies, apart from the N2 pa¬ 
rameter and the N/0 abundance ratio. N2 is clearly larger in 
barred late-type (by -0.12 dex) and lower bulge-mass (by -0.13 
dex) galaxies compared with unbarred galaxies, while this differ¬ 
ence is smaller (-0.07 dex), but still signihcant, in the early-type 
and higher bulge mass sub-sample. These differences between 
barred and unbarred galaxies in N2 translate into a larger N/0 
abundance ratio in barred galaxies by -0.12 dex in low bulge- 
mass galaxies, and -0.09 dex, in early-type, late-type and higher 
bulge-mass sub-samples. 

Finally, if we restrict our comparative analysis of barred and 
unbarred galaxies to pure star-forming galaxies, we hnd the same 


http ://data. galaxy zoo .org 
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±.U Z.U Z.J U.J ±.U Z.U 

[SII]A6717/[SII]A6731 [SII]A6717/[SII]A6731 log Up^u log Up^u 



Fig. 13. Comparative histograms of (from top to bottom and from left to right): the Balmer extinction at the emission line, c{liJ3), N2, log 
R 23 , decimal logarithm of the SFR per unit area, [S ii]/l6717/[S ii]/l6731 line ratio, logarithm of the ionisation parameter, N/0 abundance ratio and 
oxygen abundance for barred (yellow, hatched) and unbarred (purple) galaxies separately for early- and late-type galaxies. The separation between 
early- and late-type galaxies was carried out using the B/D luminosity ratio as described in Sect.j^for the non-AGN sub-sample of galaxies. 


result as with the non-AGN galaxies, i.e. late-type galaxies dif¬ 
fer in central properties, depending on whether they have a bar 
or not, with barred galaxies exhibiting larger I^sfr, ciHjS), N2, 
log (N/0) and log U, and lower [S ii]/16717/[S ii]/16731 than un¬ 
barred galaxies. Also, similar to the non-AGN sample, the dis¬ 
tribution of 12 + log (0/H) is different in barred and unbarred 
late-type galaxies and those with smaller bulge mass, with a 
marginally larger oxygen abundance in unbarred compared to 
barred galaxies. The results for late-type galaxies in both the 


pure star-forming and non-AGN samples, only differ for log R 23 . 
This shows a different distribution for barred and unbarred galax¬ 
ies in the SF sample (although with identical median values). 
However, in the non-AGN sample, the distributions are identi¬ 
cal in barred and unbarred galaxies (R-value 58%). The num¬ 
ber statistics are poorer for SF galaxies, but this reveals that the 
inclusion of the transition objects does not determine our main 
results. We will comment further on this in Sect. [121 
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Table 3. Statistics for the sub-samples of early- and late-type barred and unbarred galaxies. 




Early-type 
(based on B/D) 



Late-type 
(based on B/D) 



P-value 

barred 

median 

N 

unbarred 

median 

N 

P-value 

barred 

median 

N 

unbarred 

median 

N 

c(m 

11 . 1 % 

0.55+0.08 

26 

0.47+0.08 

43 

1.5% 

0.49+0.04 

64 

0.43+0.04 

82 

log R 23 

89.0% 

0.40+0.06 

28 

0.38+0.06 

44 

57.9% 

0.26+0.04 

64 

0.25+0.05 

73 

N2 

1 . 1 % 

-0.24+0.03 

28 

-0.31+0.04 

48 

0 . 001 % 

-0.29+0.03 

67 

-0.41+0.01 

87 

log '^SFR 

6.7% 

- 0 . 8 + 0.2 

30 

-1.3+0.2 

50 

0.92% 

-1.17+0.12 

67 

-1.33+0.12 

88 

[Sii]i6717/[S ii]^6731 

2 . 0 % 

1.25+0.04 

25 

1.31+0.05 

40 

3.7% 

1.27+0.04 

64 

1.31+0.04 

74 

logU 

85% 

-3.08+0.05 

28 

-3.08+0.03 

50 

0 . 8 % 

-3.08+0.03 

64 

-3.17+0.02 

85 

log (N/O) 

8.9% 

-0.50+0.05 

26 

-0.59+0.05 

46 

0.05% 

-0.49+0.04 

64 

-0.58+0.04 

82 

12+log (0/H) 

76% 

8.67+0.02 

28 

8 . 68 + 0.02 

50 

0.7% 

8 . 68 + 0.01 

64 

8.70+0.01 

85 


Notes. Median values and errors (95% confidence interval of the median) of the distributions of c{]Ap), R 23 and N2 gas metallicity indicators, star 
formation rate per unit area, [S ii]/l6717/[S ii]/16731 line ratio, logarithm of the ionisation parameter, N/0 abundance ratio and oxygen abundance 
in the centres of barred and unbarred galaxies in the sub-samples of early and late-type galaxies as defined in Sect. and number of galaxies in 
each sub-sample. Only non-AGN galaxies are considered. The P-values of the ^-sample A-D test for the comparison of distributions for barred 
and unbarred galaxies is shown for all cases. 

Table 4. Same as Tables and [^except for the sub-samples of low and high stellar mass bulges. 




^^bulge ^ 

10 ^-^ 

Mo 



^^bulge ^ 

10 ^-^ 

Mo 



P-value 

barred 

median 

N 

unbarred 

median 

N 

P-value 

barred 

median 

N 

unbarred 

median 

N 

c(m 

0 . 2 % 

0.46+0.03 

50 

0.42+0.04 

68 

58% 

0.52+0.07 

54 

0.50+0.06 

72 

log R 23 

49% 

0.31+0.05 

49 

0.25+0.05 

56 

67% 

0.34+0.06 

56 

0.35+0.05 

74 

N2 

0 . 001 % 

-0.30+0.03 

54 

-0.43+0.02 

68 

3.0% 

-0.26+0.03 

56 

-0.33+0.04 

84 

log 'Esfr 

0.45% 

-1.15+0.13 

54 

-1.37+0.12 

68 

2 . 6 % 

-0.99+0.17 

58 

-1.29+0.16 

87 

[Sii]i6717/[S ii]^6731 

0.97% 

1.26+0.03 

53 

1.34+0.05 

57 

2 . 0 % 

1.25+0.05 

51 

1.31+0.03 

70 

logU 

2.9% 

-3.08+0.03 

51 

-3.18+0.02 

65 

54% 

-3.07+0.03 

56 

-3.09+0.03 

86 

log (N/O) 

0.3% 

-0.51+0.04 

52 

-0.63+0.04 

63 

0.3% 

-0.47+0.04 

53 

-0.56+0.03 

80 

12+log (0/H) 

0.9% 

8.69+0.01 

51 

8.70+0.02 

65 

43% 

8.67+0.02 

56 

8.67+0.01 

86 


In summary, we observe statistically significant differences 
in the distributions of parameters for barred and unbarred galax¬ 
ies in our later-type galaxy sub-sample, regardless of whether we 
select them according to their predicted T-type (from their B/D 
light ratio) or their bulge mass. These differences appear in all 
parameters except in the metallicity tracer R 23 . Median values 
are only marginally different in most parameters (except in N2 
and log (N/0)), but they point towards a larger pR, c(I^), N2, 
log (N/0) and log U, and lower [S ii]T6717/[S ii]T6731 in barred 
galaxies. Results for earlier-type galaxies depend on whether we 
select them according to B/D light ratio or bulge mass. How¬ 
ever, with both criteria for dividing, the distributions of N2 and 
[S ii]/I6717/[S ii]/l6731 are significantly different between barred 
and unbarred galaxies. The difference in the median value of N2 
for barred and unbarred galaxies is slightly smaller than in later- 
type galaxies but is still significant. 


10. Barred versus unbarred difference trends with 
buige and totai steiiar mass 

Previous authors find observational evidence that massive barred 
galaxies have a higher current central-star-formation rate than 
unbarred galaxies of the same stellar mass ( [Ellison et al.|201l| . 
However, we see above that the effect of bars in galaxy centres 
seems to be stronger or more visible in galaxies with lower bulge 
mass, even when the barred and unbarred sub-samples have been 
created to have indistinguishable total stellar mass distributions. 


To better understand our own and previous authors’ results, 
we created a series of boxplots that show the variation of all 
parameters as a function of bulge and total stellar mass. Box- 
plots are a useful representation in our case, as sample sizes 
are not big enough to explore parameter dependences using his¬ 
tograms. Figure shows boxplots for c(HyS), N2, log R 23 , fr, 
[S ii]/l6717/[S ii]/I6731, log U, N/0, and oxygen abundance as a 
function of total galaxy stellar mass (see caption for a description 
on basic boxplots features). It is clear from these plots that, for 
the most massive galaxies, there is only a significant difference 
between barred and unbarred galaxies for the N/0 abundance ra¬ 
tio. The most massive galaxies are also the galaxies with more 
massive bulges, and this result simply confirms the results from 
Sect. For less massive galaxies (M* < M©), N2 and 

'^SFR are clearly higher for barred than for unbarred galaxies, 
with non-overlapping box-notches between barred and unbarred 
galaxies in the two lowest mass intervals. Our results are, there¬ 
fore, only partially in agreement with jEllison et al.| ( |20Il| ). We 
do find a larger fibre SFR in barred galaxies above 10^^ M©, but 
we do not see differences in SFR between barred and unbarred 
galaxies above M©. 


The differences between barred and unbarred galaxies are 
undoubtedly better correlated with bulge stellar mass (see 
Fig.[T^. All properties (except R 23 , log U and 12-rlog (0/H)) 
are significantly different between barred and unbarred galaxies 
for bulge masses < M©, but differences normally increase 

towards lower bulge masses, except for in the lowest mass 
bin, where the median values are similar within errors for both 
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Fig. 15. Boxplot showing median values of bulge stellar mass for three 
different bins in total galaxy stellar mass. 


types of galaxies. The observed differences always suggest that 
barred galaxies present larger dust extinction, star formation rate 
per unit area, [Nii]/H(r line ratio, lower [S ii]/l6717/[S ii]/16731 
(indicating a higher electron density), and larger log (N/0) and 
log U in their centres. 

We also note that the SFR per unit area is also larger on av¬ 
erage for barred galaxies in the highest bulge mass bin. How¬ 
ever, as mentioned above, we do not see differences in fr be¬ 
tween barred and unbarred galaxies in the most massive galax¬ 
ies (Fig. [T^. This is probably because bulge masses of barred 
galaxies are lower than in unbarred galaxies f or the most massive 
galaxies (Fig. [13, as already pointed out by [Coelho & Gadotti| 
( 2QTT] ). 

More massive bulges are typically classical bulges, with Ser- 
sic indexes higher than ~2.5, in contrast with the so called 
pseudo-bulges, whose surface brightness profile is better fitted 
with lower Sersic indexes and which are typically less massive 
than classical bulges ( |Gadotti|2009| ). We have also analysed the 
barred and unbarred galaxy central properties as a function of 
Sersic index. Barred and unbarred galaxies with large values of 
the Sersic index (n>3.1) do not differ in their central ionised gas 
properties, apart from the SFR, which seems to be enhanced in 
barred galaxies, in agreement with the result found for barred 
galaxies with the most massive bulges. However the global pa¬ 
rameter dependences with the Sersic index and the differences 
between barred and unbarred galaxies are not as clear as with 
bulge mass or even total stellar mass. 


11. Dependence on bar parameters 

Previous sections show clear differences in the central ionised 
gas properties between barred and unbarred galaxies. These dif¬ 
ferences imply that the gas flows induced by bars are able to 
modify gas properties, especially in lower mass or lower bulge 
mass galaxies. Simulations predict a dependency between the 
gas inflow rate and the bar strength, with stronger bars supplying 
gas towards the centre at a higher rate ( |Athanassoula| 1 992} |Re | 
|gan & Teuben|2004| ). Thus, it is natural to wonder whether these 
observational enhancements are related to dynamical or struc¬ 
tural bar parameters. 

Morphological decomposition of the galaxy sample was per¬ 
formed by |GadottT| ( |2009| ) and a number of structural parameters 
for discs, bulges and bars (in the case of barred galaxies) are 
thus available. This allows for a comparison between parameters 
such as bar ellipticity (as a proxy for bar strength) and bar length 
with central ionised gas properties. Our results are presented 
and compared with previous work below. The SFR per unit area 
(I^sfr) shows a weak tendency to increase with the effective ra¬ 
dius of the bar, the bar semi-major axis normalised to the galaxy 
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disc scale length, and the bar ellipticity. However, the dispersion 
in these relations is high and the Spearman rank correlation co¬ 
efficients (ps) are 0.28, 0.20, and 0.12, respectively. For other 
parameters, such as t he bar’s Sersic i ndex (which might be re¬ 
lated to the bar’s age; |Kim et al.|2015] ) ps is lower than 0.1. This 
agrees with previous works on central SFRs from infrared emis¬ 
sion (e.g. [Pompea & Rieke||l990t [Roussel et al.||200f] ), which 
did not find evidence that stronger bars could produce a larger 
enhancement in SFRs. [Eflison et ak] ( |201 Ij ) in turn, did not find 
a correlation between bar length and fibre SFR in the central 
region of galaxies from SDSS optical spectra. However, other 
authors have reported a positive relation between central star for¬ 
mation and bar parameters, but their results are, in some cases, 
conflicting: [Oh et aL] ( |2Q12| ) report a higher effect of bars on cen¬ 
tral SF but only for the reddest galaxies where the effects are 
more pronounced with increasing bar length. In contrast, other 
authors ( jWang et~ar]|2Q12| ) And that only strong bars can en¬ 
hance central SF, but the degree of enhancement depends solely 
on bar ellipticity, and not on bar size or bar mass. This matches 
with [Zhou et al.| ( |2()15| ), who find a weak positive trend between 
SFR in bulges and an ellipticity-based bar strength parameter. 

Several studies have also tried to confirm observationally the 
effect of bars on the central oxygen abundance and on the ra¬ 
dial oxygen abundance gradient (e.g. [Florido et al.||2012| and 
references therein), but results are also conflicting. For example, 
there is disagreement on whether bars produce a larger oxygen 
abundance in galaxy centres (e.g. [Fllison et al.||2Qll|), equal to 
that found in unbarred galaxies (e.g. jChapelon et al.||1999t jCa- 


jcho et al.|2Q14| ) or whether 0/H is lower in barred galaxies than 
in galaxies without a stellar bar ( jConsidere et al.|2()00[ jPutil &| 
|Roy|199^ . However there seems to be a consensus that neither 
central nor the 0/H abundance gradient depends on bar strength 


(jChapelon et al.|199^ Considere et al.|2000 jCacho et al.|2014[ 


[Fllison et al.||201f] )r Our results agree with previous studies, 

since we also find no clear trend between 12-hlog(0/H) and bar 
parameters (effective radius, normalised length, ellipticity) with 
correlation coefficients ps lower than 0.2, with a slightly better 
correlation with the bar boxiness (ps =032). 

We also find no correlation between the N/0 abundance ra¬ 
tio and any bar parameter, with ps lower than 0.1 in all cases, 
with the exception of a weak trend for the N/0 to increase 
with the bar effective radius (ps =0.15). The c(Uj3) and the 
[S ii]/16717/[S ii]/I6731 emission-line ratio do not show any clear 
dependence on bar parameters either (ps ^0.2). 


12. Discussion 

In this paper we perform a comparative analysis of the properties 
of the ionised gas in the centres of disc galaxies that possess a 
stellar bar, with those without a bar structure. Our analysis is 
based on SDSS spectra. 

The SDSS fibre size (3") and our face-on galaxy sample 
redshift range (0.02 < z < 0.07) implies that, in this work, we 
are analysing emission of the ionised gas located within galac- 
tocentric radii of between approximately 0.6 and 2.1 kpc. The 
ionised gas-emitting regions located at these galactocentric dis¬ 
tances are usually termed as "nuclear H ii regions" or "H ii region 
nucleus" when they are located in the immediate neighbourhood 
of a galactic centre, and "hotspots" or "circumnuclear Hii re¬ 
gions" when these regions surround the galaxy centre, frequently 
arranged in a ring or pseudo-ring shape (see e.g. jKennicutt et al.] 
[1989[[Ho et al.|1997a|[Diaz et al.|2007[ ). The spectroscopic stud¬ 
ies of nuclear and circumnuclear spatially resolved H ii regions 
are still limited to a small number of targets. These have been 
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the focus of previous in-depth spectroscopic studies, mainly with 
the aim of comparing the physical and stellar content properties 
to those of normal disc H ii regions (e.g. Kennicutt et al.||1989 


Ho et al.||1997at [Bresolin & Kennicutt||2002| ), studying differ¬ 
ences in properties with the Hubble-type of the host galaxy ([Ho| 


et al. 1997a) and/or deriving chemical abundances. These nu¬ 


clear and circumnuclear Hii regions are expected to be those 
with the highest metal content because of their location in the 
inner parts of spiral discs ( |Diaz et al.|2007l|Dors et allMTT] ). To 
our knowledge no published work so far exists that compares the 
properties of resolved nuclear and circumnuclear H ii regions in 
barred and unbarred galaxies for a large number of galaxies in a 
systematic way. 

Our derived values for internal B aimer extinction, 
[Sii]/16717/[Sii]/16731 ratio, oxygen abundance, and N/0 
ratio are within the values quoted by authors analysing indi¬ 
vidual, resolved nuclear, and circumnuclear Hii regions (e.g. 


Kennicutt et al.[[1989l 

Bresolin & Kennicutt[[2002[ [Diaz et al.[ 

2007| Dors et al.||2008 

). Our finding, that the current SFR and 


electron density (as traced by [S ii]/16717/[S ii]/16731) are also 
larger i n barred systems, is also in agreement with |Ho et al.] 
( |1997b| ), who made a spectroscopic survey of the nuclei (using 
a 2"x4"aperture) of nearly 500 nearby galaxies. In particular, 
they only found a larger SFR in barred galaxies earlier than 
Sbc, compared with unbarred counterparts, with no difference 
between barred and unbarred late-type spirals. We also find a 
larger average for SFR per unit area in barred galaxies in our 
early-type sub-sample (see Fig. [T3 but the distributions for 
barred and unbarred galaxies are only slightly different, while 
the difference in both distributions is statistically different for 
the late-type galaxies, according to the A-D test. We note that 
the most significant differences between barred and unbarred 
galaxies (in all gas properties) are obtained for lower bulge or 
total stellar mass galaxies (i.e. later-type galaxies). The total 
stellar galaxy or bulge mass distribution for the galaxy sample 
of |Ho et ar] ( |1997b| ) is not shown and we cannot go any further 
in comparing both results. 


12.1. Enhanced [N ii]A6583/Ha in centres of barred galaxies 


The most striking difference between barred and un¬ 
barred galaxies obtained in this work comes from the 
observed [Nii]/16583/HQf emission-line ratio. An enhanced 
[Nil]46583/Ha' emission-line ratio in the galaxy centres was 
first found by Stauffer ( 1982|) and later by different authors (see, 
for example ^ ^ ^ ' 

[erari[20T5] ). 


■Cennicutt et ai.||1989t |Ho et al.||1997at |Sanchez| 
io et al.| ( |1997a| ) analysed the integrated spectra 


of a sample of 206 galaxy nuclei. They reported that early-type 
galaxies had the largest [Nii]/I6583/Haf. According to these 
authors, it might be an indication that nitrogen is selectively 
enriched in the centres of bulge-dominated galaxies. This agrees 
with the general trend found in this work (without distinction 
between barred and unbarred systems), since early-type galaxies 
are also, on average, more massive, and they host bulges with 
a greater mass than late-type galaxies (see Figs. [T^ and [T^. In 
this work we go a step further and show that this line ratio is 
larger in the centres of barred galaxies and that the bar effect, in 
enhancing this line ratio, is more important in galaxies with less 
massive bulges. 


A number of effects can enhance the [Nii]/I6583 to Ha 
emission-line ratio, namely the nitrogen excitation by a source 
other than photoionisation (i.e. AGN or shocks) that enhances 
the [N ii] emission, compared with emission from recombination 
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Fig. 18. Comparison of shock- and AGN-sensitive emission-line ratios 
[N ii]/l6583/HQf and [O i]/l6300/HQf for barred and unbarred galaxies. 
Galaxies that would be considered as AGN accordi ng to the [O iii]/I^ 

[O fi/Hof or [O [S ifi/Hcr diagrams and the |Kewley et al.| ( [2506] ) 

separation line are marked with blue and green circles, respectively. 


lines, and a larger abundance of nitrogen atoms compared with 
oxygen. 

We n ow analyse these cas es separately. We have used the 
criteria of |Kewley et al.| ( [2Q0T] ), based on the [O iii]/I5007/HyS vs. 
[Nii]/I6583/Hcr, to remove AGN from the galaxy sample. This 
separation criteria allows for objects of the composite area to be 
included, which may add a number of galaxies in which shocks 
or an AGN can enhance the [N ii]/I6583 emission, but also pure 
photoionised objects with high N/0 ratio ( [Perez-Montero & 


|Contini|[2009| ). We note that pure SF galaxies also show a sig¬ 
nificant difference in [Nii]/I6583/Haf between barred and un¬ 
barred ga]uxies_(s^Tab]^J^. However, to check whether the 
use ofr 


Kewley et al. 


(2001|) criteria includes a high number of 


shock-ionised barred galaxies, and enhances the observed dif¬ 
ferences (e.g. in the right-hand panel of Fig. between barred 
and unbarred galaxies, we created Fig. This figure shows 
[Nii]/I6583/Haf as a function of the shock and AGN sensitive 
[Oi]/16300/Ha' line ratio. We can see that galaxies with large 
log([Oi]/16300/Ha') in their centres (> -1) are also galaxies with 
high [Nii]/I6583/Haf and, according to the [O iii]/I5007/H/5 vs 
[Qi]/163 00/Ha' BPT diagram , those galaxies would be classified 
as AGN ( [Kewley et al.|2006] circled in blue). However, there is 
roughly the same number of barred (11) and unbarred (10) galax¬ 
ies in this area of the plot. Of those galaxies, only three barred 
galaxies and two unbarred galaxies would have been classified 
as AGN with the [Oiii]/15007/HyS vs [ S n]/l/16717,673 l/Ho f dia¬ 
gram and the separation line given by [Kewley et al.| ( |2006| ), cir¬ 
cled in green. The plot clearly shows that even by removing those 
galaxies, the average [Nii]/16583/Haf is larger in barred than in 
unbarred galaxies. 

In addition, if shocks or AGN excitation were responsi¬ 
ble for a larger [Nii]/16583/Haf in barred galaxies, we would 
also expect to see, on average, larger values in other AGN 
and shock-sensitive emission-line ratios, such as [O i]/16300/Haf 
or [Sii]/l/16717,6731/Haf, in barred galaxies when compared 
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to their unbarred counterparts. The distributions of the three 
emission-line ratios are shown in Fig. It is clear that 
the distributions are only significantly different in the case 
of [Nii]/16583/H(r, while, in the case of [Oi]/I6300/HQf and 
[S ii]/I/16717,673 l/Hof, the mean values for barred and unbarred 
galaxies are identical, and that the distributions are not statisti¬ 
cally different, with A-D test P-values greater than 30%, com¬ 
pared to 0.0009% for [Nii]T6583/Ha. 

A difference in N/0 abundance ratio in the centres of barred 
and unbarred galaxies is thus the most probable explanation for 
the observed enhancement in the [Nii]/I6583/H(r emission-line 
ratio in barred galaxies. But what is responsible for an over¬ 
abundance of nitrogen compared with oxygen in barred galax¬ 
ies? 

The N/0 abundance ratio is a useful parameter to understand 
chemical evolution of galaxies (see e.g . [Vila Costas & Edmunds 


1993[|Perez-Montero & Contini|2009| ), as both elements are pro- 


to the galaxy centre will depend on the bar strength (e.g |Mar- 
|tin & Friedli||1999} [Regan & Teuben|[2004| ). The gas flow in 
the bar gravitational potential will form shocks along the bar’s 
leading edges, givin g way to gas flows towards the centre (e.g. 
|Athanassoula|1992] ). The availability of gas in the disc will also 
alter the amount of gas accreted into the bulge region and, as the 
bar evolves and the gas reservoir diminishes, the inflow rate will 
also decrease (jAthanassoula et al.| |2005| ). Therefore strong gas 


duced by different mechanisms and by stars of different masses. 
While oxygen is always primary and created in massive stars, 
nitrogen is synthesised via the CNO and CN cycles in H burning 
stars, and is mostly secondary (the seed C and O were already 
present in the placental cloud when the star formed) for metal- 
licity higher than 12-l-log (0/H) = 8.3, while for low metallicity 
very little secondary N is produced. Also, N can be produced in 
stars of all masse s, but mainly in int ermediate mass (between 4 
and 8 M©) stars ( [Henry et al.||2000| ). Consequently, in a single 
burst of star formation, we expect that most of the N will be re¬ 
leased to the ISM after O, with a time-delay roughly equal to the 
stellar lifetime of the main producers of N (intermediate mass 
stars), i.e. from -100 to 500 Myr after the oxygen was released. 
Star formation in spiral galaxies m ight be a continuous process 
(e.g. [Gavazzi & Scodeggio[[1996[ ) which complicates this sce¬ 
nario. In addition, gas flows (within the galaxy and from outside) 
can also, in theory, alter the initial abundance ratios. 

Both chemical evolu tion models ([Molla & Gavilan[[2010[ ) 
and observational works ( [Mallery et al.|2007[ ) show that the N/0 
abundance ratio depends on the star formation history. In simple 
terms this means that, in a closed-box model and at a fixed metal¬ 
licity, all that is required for a galaxy or region to increase its N/0 
is for its current SFR to be less than its past average SFR, so that 
comparatively fewer high mass stars are now formed. Thus, al¬ 
though newly synthesised oxygen is being released into the ISM 
from any massive stars that are currently forming, intermediate- 
mass stars from previous generations (when SFR was higher) 
have had time to evolve and a larger amount of N has been re¬ 
leased into the ISM, increasing the N/0 abundance ratio. How¬ 
ever, metal-poor gas fiows have a relatively small effect on the 
N/0 abundance ratio according to some models (e.g. [Edmunds] 
[1990[ ). While the oxygen abundance is significantly decreased 
during infall, as a result of dilution of the present gas, the N/0 
ratio is much less sensitive. The mass infall must be much larger 
than the galaxy mass (or bulge mass in our case) and the infall 
rate must be higher than the SFR to have any large deviations 
from the expected closed-box behaviour ([Koppen & Hensler[ 
[2005] ). 

A different history of SF in the centre of barred galaxies 
compared with unbarred galaxies is expected from numerical 
simulations and can, at least, qualitatively explain the higher 
N/0 abundance ratio observed in bulges less massive than - 
10^-^ - 10^^ Mq. Estimating the efficiency of the bar in trans¬ 
porting material to the centre, using numerical simulations, is 
not easy because the inflow rate strongly depends on the simula¬ 
tion and the chosen parameters (e.g. [Kim et al.|20I5| . However, 
it is expected that the efficiency of the bar, in driving material 


fiows can be expected soon after the bar forms, which presum¬ 
ably leads to SF in the inner parts (e.g. [Cole et al.|20I4[ ). Later 
on, the bar weakens and gas fiows continue to transport mate¬ 
rial towards the central regions, where it can form stars, but at 
a lower rate. Thus, according to chemical evolution models, this 
could explain the high N/0 in barred galaxies compared with un¬ 
barred galaxies. Apart from bars, other mechanisms can produce 
gas fiows towards the centre of the galaxies (holding or not a 
bar), for example tidal forces by an interacting companion (e.g. 
Cox et al. 2008) or asymmetr ic structures like spiral arms (e.g. 
[Kormendy & Kennicutt[2004} , but these might be more sporadic, 
and not produce a continuous flow, as in the case of barred galax¬ 
ies. 

Given that a considerably strong bar is present, late-type 
galaxies have a large reservoir of gas and can, therefore, sup¬ 
ply a large amount of gas towards the centre. This is not the 
case for bulge-dominated galaxies and this might be, at least 
partially, one of the reasons why differences between barred and 
unbarred galaxies are smaller or non-existent in galaxies with 
massive bulges (> 10^^ M©). The fact that we do not find differ¬ 
ences in central N/0 abundance ratio in the most massive bulges 
could indicate that bar effects on central gas properties are either 
negligible (maybe because there is not much gas in the disc to 
feed the inner parts) or invisible. This may be because there is 
another, more important and dominant, mechanism that masks 
the expected bar effects. 

In any case, for massive bulges (which are also the ones with 
larger N/0 abundance ratio, see Fig. [nj it would be more com¬ 
plicated to detect an increase in N/0 in barred galaxies. The rea¬ 
son being that, given this time-delay scenario in the ejection of 
N, at high abundances the increase in the N/0 abundance ra¬ 
tio is smaller ( [Coziol et al.[[I999T ). However, the fact that we 
do not see differences in other properties (c(HyS), 12-rlog (0/H), 
log U), or that these differences are smaller (as with fr oi* 
[S ii]/I6717/[S ii]/I6731) in the high bulge mass range, points to¬ 
wards a different bulge origin, formation, or evolutionary pro¬ 
cess, that is dependent on bulge mass (or total galaxy stellar 
mass). 

Bulges are usually classified as two different types: disc-like 
bulges or pseudo-bulges, and classical bulges. The former are 
generally less massive than classical bulges and are more com¬ 
mon among late-type galaxies. The formation of classical bulges 
is thought to be faster and normally externally driven, by means 
of gravitation al collapse or hi erarchical merging o f smaller ob¬ 
jects (see e.g. [Gadotti|2009[ or [Athanassoula|200>5] for details on 
the formation mechanisms of the different types), while disc-like 
bulges are thought to be formed by internal secular processes, 
predominantly driven by bars (e.g. [Kormendy & Kennicutt[2004t 
[Athanassoula|2005] ). 

There is observational evide nce for a delay in the bar for¬ 
mation for less massive galaxies ( [Sheth et al.|2QI2[ ). Therefore, 
both bars and bulges in massive galaxies, as mentioned above, 
form at earlier times. These bulges might be classical as they 
form in more violent environments, with more gas and larger 
merger rates. Moreover, recent numerical simulations ( [Kraljic 
et al.|20I2| ) show that the formation of bulges in spiral galaxies 
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is influenced by the bar presence. They find an over abundance of 
star formation in the bulge for barred galaxies after the bar forms. 
It is thus possible that the difference observed in the central stel¬ 
lar metallicities of massive galaxies ( [Coelho & Gadotti]|2011| ) 
is set at the time when the bar forms, and consumes and redis¬ 
tributes most of the gas. This would be compatible with the fre¬ 
quent co-existence of both type of bulges in the centres of galax- 
ies as observed in previous works ([Peletier et al.|20Q7t [Mendez- 1 

[Abreu et al.|2014[ [Erwin et al.|2015| ). In this scenario, low mass 

systems do not develop a massive classical bulge, and the bar 
will form at later times. Furthermore, the central gas accretion 
lasts longer, secularly forming disc-like bulges in barred galax¬ 
ies. This would generate a difference in the ionised gas properties 
between barred and unbarred low-mass systems as observed in 
the present work. 

Therefore, our result, that central properties of ionised gas 
vary predominantly in galaxies with lower mass bulges, gives 
strength to a different origin or evolutionary process for less and 
more massive bulges, with the gaseous phase of the former cur¬ 
rently having a closer relation with bars. 

12.2. Comparison with recent works 

Several works have tried to find observational evidence for the 
effect of gas flows by analysing the gas content, SFRs and/or 
metallicities in galaxy centres (e.g. [Ho et al.|[1997bt [Martin &] 
[Friedli[p~999[ ) but their results were inconclusive. This was, in 
part, because these works were carried out with small data sam¬ 
ples and sample selections may have introduced biased results 
caused by, now well-known, relations between the analysed gas 
properties with Hubble type and/or total stellar mass. In an effort 
to s ide-step these proble ms, more recently [Ellison et ar] ( [201I[ ) 
and [Cacho et~ar] ( [2014[ ) have comparatively analysed gaseous- 
phase metallicities and SFRs in the centres of barred and un¬ 
barred galaxies in large samples (from -900 to 1500 galaxies, 
respectively) by using SDSS spectra. In both cases the galaxy 
sample s have been taken from the catalogue of [Nair & Abraham] 


201 Op requiring axial ratios bja > 0.4 (i.e. i < 68 °, compared 


mann et ak] ( [2003 a[ ) to separate star-forming dominated galax¬ 
ies from active galaxies, which rejects galaxies from the tran¬ 
sition area, where pure star-forming galaxies with high nitro- 
gen abundance are expected to be lo cated (e.g. [Perez-Montero 


& Contini|2009[). Ellison et al.[ ( [2011| ) and the present study use 


Kewley et al.[ p001[ ) and, therefore, include transition objects 


to 26° in our sample) and redshift range z < 0.06 and z < 0.1 in 
Ellison et~ar] ( [201l| ) and [Cacho et al. ( 2014[ ), respectively. Spe- 


Bars are visually classified in the Nair & Abraham (2010) catalogue 


in our analysis. Our samples are thus e xpected to contain m ore 
nitrogen-rich galaxy centres. In addition, Ellison et al. ( 2011|)use 
the [Nii]/[Oii] empirical calibration ( [Kewley & Dopita[[2002[ ) 
to obtain 12-rlog (0/H), which depends on the N/0 abundance. 
We believe their reported larger 12-rlog (0/H) in barred galaxies 
might indicate a larger N/0 ab undance ratio in bar red galaxies 
compared to unbarred galaxies. [Cacho et ar] ( [2014| ) use the R 23 
parameter to obtain oxygen abundances, but later on, in an effort 
to explain the differences wit h Ellison’s results, th ey apply the 
same empirical calibration as [Ellison et al.[ ( [ 2 QTT] ), [Nii]/[Oii]. 
With this calibration they also find larger I2-rlog (0/H) in barred 
galaxies compared t o unbarred galaxies, b ut the difference is not 
as significant as in [Ellison et al.[ pQll[ ) because their sample 
lacks a significant number of nitrogen-rich galaxies, as explained 
above. We also get larger 12-rlog (0/H) in barred galaxies when 
using empirical calibrations involving [Nii]/I6584, such as the 
on es given by[Pettini & Pagel] ( 2004[ ) ; [Kewley & Dopita[ ([20Q2[) 
or [Bresolin ( 2007[ ), but no difference when R 23 or the [Perez- [ 
[Montero[ ( |20i4[ ) method is used, as described in Se ct.[ 6 j 


Finally, we would like to mention the fact that Ellison et al 


( [ 201 1 | ) only find differences in the most massive galaxies of their 
sample (M* > 10 ^^ Mq) does not contradict our results, since 
our galaxy sample contains only galaxies more massive than 
10^^ Mq. Our results do not match in that we do not find differ¬ 
ences between barr ed and unbarred gala xies more massive than 


Mn, while 


Ellison et al. 


cial care was taken by the authors to avoid biases due to dif¬ 
ferences in inclination, redshift and stellar mass between their 
barred and unbarred galaxy samples. However, and in spite of 
the similarity in their sam ples and in the data an alysis, they arrive 
at differing results. While [EITison et al.| ( [ 2 Qll| ) find that the SFR 
and the oxygen abundance in the galaxy centres are considerably 
larger in barred galaxies more massive t han 10 ^^ Mq compa red 
with unbarred galaxies of the same mass, [Cacho et al.[ ( |2014[ ) do 
not find statistically significant differences between barred and 
unbarred galaxies, only a trend (which decreases with increas¬ 
ing mass) in their early-type sub-sample, for barred galaxies to 
exhibit a larger 12-rlog (0/H) than unbarred galaxies. 

The study presented h ere is in apparent disagr eement with 
both [Ellison et al.[ ( |2011[ ) and [Cacho et ak] ( [2014| ), as we have 
found significant differences between barred and unbarred galax¬ 
ies, but we did not see differences in the oxygen abundance be¬ 
tween barred and unbarred galaxies. However, our results can ex¬ 
plain, to a large extent, the disagreement of [Ellison et al.[ ( [MTT[ ) 
and [Cacho et al.[ ( |2014[ ) results, both with each others’ and with 
ours. This mostly arises from the way in which star-forming 
galaxies are selected, and from the empirical calibration used 
to obtain oxygen abundances. [Cacho et al.[ ( [2014[ ) use [Kauff-[ 


( [201 Ij ) finds larger SFRs and 
12-rlog (0/H) up to their upper limit mass (-10^^-^ M©). 

We believe that the results presented here show the foot¬ 
prints of bar-induced secular evolution of galaxies in the central 
gaseous component, which seems to be either more efficient or 
more relevant in galaxies with less massive bulges (i.e. later-type 
galaxies), at least regarding the current central properties of the 
gas. 


13. Summary and conclusions 

We analysed the properties of the ionised gas in the centre of 
a sample containing 251 and 324 barred and unbarred face- 
on galaxies, respect ively (173 and 265 after removing AGN, 
Kewley et al. 2001) of total stellar masses between 10^^ and 
10"-^ Mq, where both sub-samples have the same total galaxy 
stellar mass and redshift distributions. We compared the distri¬ 
butions of internal extinction at HjS, SFR per unit area (Zs fr), 
[Sii]/I/I6717,6731 emission-line ratio, empirical tracers of oxy¬ 
gen abundance (R 23 and N2=[Nii]/l6854/Haf), oxygen abun¬ 
dance, N/0 abundance ratio, and ionisation parameters, for 
barred and unbarred galaxies separately. Our main conclusions 
are summarised below. 


- The whole non-AGN barred galaxy sample has statis¬ 
tically significant different distributions of c(Uj3), Zsfr, 
[Sii]T/I6717,6731, N2, and log(N/0) in their centres com¬ 
pared with unbarred galaxies. The differences are remark¬ 
ably evident in N2 and log(N/0), which are ~0.1 dex larger 
in barred than in unbarred galaxies. Barred galaxies also tend 
to have larger central internal extinction, Zsfr, and electron 
density than unbarred galaxies, but in these cases the differ¬ 
ences in median values between both distributions are within 
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errors. We do not find differences between barred and un¬ 
barred galaxies, neither in the central value of the R 23 pa¬ 
rameter, nor in 12-l-log (0/H), as determined from the grids 
of photoionisation models by |Perez-Montero| (|20 1 4 } . 

- When we split our sample into earlier and later types (ac¬ 
cording to their B/D light ratio or their bulge mass), the 
differences increase towards later-type galaxies. Especially 
significant is the difference in N2 that is ~0.13 dex larger 
in late-type barred galaxies compared to unbarred galaxies, 
against a difference of ~0.07 dex between barred and un¬ 
barred galaxies in the earlier-type sub-sample. This differ¬ 
ence in N2 translates into a difference in N/0 abundance ra¬ 
tio of 0.09 - 0.12 dex from photoionisation models. 

With this separation into earlier and later types, the oxygen 
abundance distribution appears to differ in barred and un¬ 
barred late-type galaxies, with a slightly larger median value 
for unbarred galaxies. There is no difference in 12-l-log (0/H) 
between barred and unbarred galaxies in the early-type sub¬ 
sample. In addition to 12-rlog (0/H), barred and unbarred 
early-type galaxies seem to have similar central c(Uj3), R 23 
and log U. The distributions of fr and N/0 are different 
in barred and unbarred early-type galaxies (indicating larger 
median values in barred galaxies) depending on the way we 
separate early- from late-type galaxies. 

- The parameter differences between barred and unbarred 
galaxies increase for galaxies with lower mass systems and 
seem to correlate better with bulge mass, except for the SFR 
per unit area, which seems to have a stronger dependence on 
total galaxy stellar mass. 

- We find no clear relation between the central ionised gas 
properties in barred galaxies and the bar structural param¬ 
eters (effective radius, bar length, ellipticity or bar Sersic 
index), but we do find weak trends for the central I^sfr to 
increase with bar effective length. 

- The differences observed between barred and unbarred 
galaxies can be explained qualitatively with current knowl¬ 
edge of bar evolution simulations and chemical evolution 
models. 

Taken all together, our results imply that the presence of bars 
does alter the properties of the central ionised gas of galaxies, in¬ 
creasing their electron density, internal extinction, SFR per unit 
area, and, very notably, the N/0 abundance ratio. These bar ef¬ 
fects on central ionised gas are more visible in later-type galax¬ 
ies, or more precisely, in those galaxies with bulges that are less 
massive than ~ 10^^ M©. 

This work presents some of the clearest observational evi¬ 
dence so far for the effects of bars in galaxy centres. It also em¬ 
phasises (1) the usefulness of the N/0 abundance ratio, which 
gives important clues on star formation history of galaxy cen¬ 
tres, and (2) the risk of using certain empirical calibrations for 
oxygen abundance estimations, which can mask a high N/0 un¬ 
der a false high O abundance. 
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Appendix A: Bulge-to-disc flux ratios for different 
galaxy morphological types. 

The decrease of the B/D flux ratio along the Hubble sequence 
from early- to late-type galaxies is well known, but this ratio, 
among others, depends on the structural component-fitting pro- 
cedure of the galaxy, on inclination, and on dust-extinction (see 
Graham & Worley [2008 1 and references therein). Morphological 
classification is only available from the largest catalogue so far 
( Nair & Abraham|2010| ) for less than 50% of our galaxy sam ple. 
Our aim here is to use the B/D fiux ratios ( |Gadotti||2009| ) and 
the morphological classification in T-types for our sub-sample 
of galaxies included in the [Nair & Abraham] ( |2010| ) catalogue, 
to separate our galaxies into early- and late-types, according to 
their B/D fiux ratios. 

Figure [at] shows the dependence of the logarithm of the g-, 
r-, and /-band B/D light ratios ( from BUDDA, |Gadotti||2009| ), 
with the corresponding T-type ( |Nair & Abraham] |2010| ) for 



T-type (NAIO) 


Fig. A.2. Decimal logarithm of the bulge-to-disc flux rati o in the /-band 
as a functio n of the morphological T-type, as given in the |Nair & Abra-| 
|ham| ( |2010| ) catalogue for the galaxies of our sample included in this 
catalogue. The colours r epresent different val ues of the predicted T-type 
Galaxy Zoo 2 parameter | Willett et~^ ( |2013| ). 


barred and unbarred galaxies occurring in both samples simul¬ 
taneously. The top and bottom left-hand panels show this rela¬ 
tion for the three different photometric bands, with small crosses 
for unbarred galaxies and small, filled circles for barred galax¬ 
ies. The large crosses and big empty circles indicate the median 
log (B/D) value for each T-type value for unbarred and barred 
galaxies, respectively. As can be seen, there are no clear differ¬ 
ences between barred and unbarred galaxies. The bottom right- 
hand panel shows the median log (B/D) values, for barred and 
unbarred galaxies together, for each T-type for all photometric 
bands. The differences in the median B/D light ratio between the 
different bands are, in general, smaller than the error bars. 

Our median values correspond well (within errors) to those 
from [Graham & Worley | ( [200^ their Fig. 6) in the T-type range 
in common (between T-type 0 and 5), except for T=0, where we 
obtain a marginally larger B/D ratio in the r and / bands. These 
differences cou ld arise from poor statistics at T=0 in [Graham &| 
[Worley[ (2008[ ) (only three galaxies versus 30 in our sub-sample 
of galaxies are included in the [Nair & Abraham[ ( [201Q[ ) cata¬ 
logue). 

We point out that a number of unbarred galaxies in our 
sample (~15), supposed to be disc galaxies, are classified with 
T-type<-3 by [Nair & Abraham] ( [2010| ), i.e. as ellipticals (see 
Fig. [A.l[ or [A.2[ ). These galaxies are very difficult to classify, 
since a face-on disc without spiral arms or a bar looks just like a 
round elliptical. [Gadotti[ (2009[ ) classified galaxies after looking 
at isophotal contours and the profile very carefully, but as stated 
in the cited paper, the estimated misclassification frequency is 
~5-10%. With these uncertainties in mind, we have opted to 
keep these galaxies in our sample, and to consider them as SO 
disc galaxies. 

The B/D fiux ratios decrease towards late-type spirals as ex¬ 
pected. We have used this trend to create sub-samples domi¬ 
nated by early- or l ate-type spirals. In our su b-sample of galax¬ 
ies included in the [Nair & Abraham[ ( [2Q10[ ) catalogue, a value 
of log (B/D)=-0.30 (or B/D=0.5) more or less separates galax¬ 
ies with T-type<2 (earlier than Sa) from those later or equal to 
T-type=2 (Sab), but the contamination from late-types (with the 
previous definition) in the early-type sub-sample and vice versa 
is approximately 12% and 30%, respectively. The predicted T- 
type, as given in the Galaxy Zoo 2p]catalogue (hereinafter GZ2 
[Willett et al.[[2013] ), allows us to improve this separation. We 


http ://data. galaxy zoo .org/ 
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then consider early-type galaxies as those galaxies of the sample 
with either B/D>0.5, or those having simultaneously B/D<0.5 
and T-type, as predicted by GZ2, lower than -1. Similarly, we 
consider as late-type galaxies those with B/D<0.5 and a T-type 
as predicted by GZ2, larger than or equal to one. In this way we 
are able to include ~70% and 90% of the early (T-type< 2) and 
late types (T-type >2) in the sub-samples of early- and late-type 
galaxies, respectively, and the contamination from late (early- 
types) in the early (late-type) sub-sample is -10% (25%). These 
numbers are the same in all three bands. However, we have used 
the /-band, as the B/D flux ratio is available for all sample galax¬ 
ies in this filter. 
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Fig. 14. Comparative histograms of (from top to bottom and from left to right): the Balmer extinction at the Yip emission line, R 23 , 

N2, logarithm of the ionisation parameter, logarithm of the SFR per unit area, the [S ii]/l6717/[S ii]/l6731 line ratio, logarithm of the ionisation 
parameter, N/0 abundance ratio, and oxygen abundance for barred (yellow, hatched) and unbarred (purple) galaxies separately for non-AGN 
galaxies with bulge mass lower than 10^ ^ M© and for galaxies with heavier bulges (> 10^ ^ M©). 
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Fig. 16. Boxplots showing the distribution in N2 =log([Nii]6583/HQf), R23, SFR per unit area and [S ii]/l6717/[S ii]/l6731 line ratio, 

decimal logarithm of the ionisation parameter, N/0 abundance ratio and oxygen abundance, as a function of the total galaxy stellar mass (M^) for 
barred (yellow boxes) and unbarred (purple) galaxies. Each box represents the data distribution for galaxies in a 0.42 dex interval of logarithm in 
stellar mass centred at the abscissa axis values of the box position. Inside each box, the central horizontal black straight line marks the median 
value. The lower and upper quartiles are represented by the outer edges of the boxes, i.e. the box length is equal to the inter-quartile range, and 
therefore the box encloses 50% of the data points. The notches mark the 95% confidence interval for the median value. The whiskers extend to 
the most extreme data point within 1.5 times the inter-quartile range. Galaxies that do not fall within the reach of the whiskers are considered 
as outliers (yellow and purple circles for barred and unbarred galaxies, respectively). The number of galaxies in each box, in order or increasing 
stellar mass, ranges from 35 to 39, 55 to 59, and 13 to 14 for barred galaxies, and from 57 to 64, 58 to 70, and 15 to 21, for unbarred galaxies. 
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Fig. 17. Same as Fig. Il^ but as a function of the bulge mass. Each box represents the data distribution for galaxies in a 0.85 dex interval of 
logarithm in bulge stellar mass centred at the abscissas values of the box position. The number of galaxies in each box, in order of increasing 
stellar bulge mass, ranges from 9 to 10, 66 to 71, and 28 to 32 for barred galaxies, and from 18 to 22, 70 to 81, and 28 to 32, for unbarred galaxies. 
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Fig. 19. Histograms showing the distribution of AGN and shock-sensitive emission-line ratios (from left to right: [N ii]/l6583/Hcr, 
[Sii]/l/16717,6731/HQf, and [O i]/l6300/HQf) for barred and unbarred galaxies separately. All galaxies classified as non-AGN are included (see 
Sect.j^. The barred/unbarred distributions are only different for [N ii]/16583/Hq', and therefore the source of the differences in that line ratio must 
originate from something other than shock ionisation and/or AGN contamination alone. 
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Fig. A.l. Logari thm of the B/D flux ratio in the g, r, and /-band versus morphological type (expressed as T-type) for the galaxies of the sample 
contained in the |Nair & Abraham| ( |2010| ) catalogue. Big crosses and big empty circles mark the median log (B/D) for each T-type value for 
unbarred and barred galaxies, respectively. The bottom right-hand panel shows the median log (B/D) in each band for each T-type value. Error 
bars represent the 95% confidence level of the median value for the distribution of B/D values at each T-type value. The blue solid straight line 
marks B/D=0.5, that we have used to separate early from late-type galaxies. 
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